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INTRODUCTION 

Each subroutine in the Scientific Subroutine Package is documented by means of 
a flowchart. In the guide to the flowcharts each flowchart is listed under a 
general heading: statistics, special matrix operations, matrices, or other 
mathematics; and under a subheading relating to its basic function. 

SYMBOLS USED 

Figure 1 illustrates the various blocks used in the flowcharts and their parti- 
cular meaning. Lines connecting these blocks are made up of asterisks. 
Arrows showing the direction of flow are represented by a "V" for downward, 
a left parenthesis " (" for right to left, and a right parenthesis ")" for left to 
right. Special symbols are used for mathematical operators within the flow- 
chart boxes. These are "*", "&", and "-" for multiplication, addition, and 
subtraction respectively. 

BLOCK NUMBERING 

Numbers appearing inside the connector circles refer to the chart block 
location on the page; these are numbered vertically zero to nine and horizon- 
tally A to C. Numbers shown outside the blocks correspond to FORTRAN 
statement numbers in the programs. 

EQUATIONS 

Many of the flowcharts contain references to equation numbers. These equations 
can be found in the mathematical description for the particular subroutines in 
the Programmer's Manual. 



czz> 



Enter or exit block 



Processor block 



Modification block 




Decision block 




Call to subroutine block 



Connector to block within page 



Connector to block on another page 



Figure 1 . Flowchart blocks 



GUIDE TO FLOWCHARTS 
FLOWCHARTS 

Statistics 

Data Screening 

TALLY — totals, means, standard deviations, 

minimums, and maximums 8 

BOUND — selection of observations within bounds 9 

SUBST — subset selection from observation matrix 10 

ABSNT — detection of missing data. 11 

TAB1 — tabulation of data (1 variable) 12 

TAB2 — tabulation of data (2 variables) 13 

SUBMX — build subset matrix 15 

Elementary Statistics 

MOMEN — first four moments 16 

TTSTT — tests on population means 17 

Correlation 

CORRE — means, standard deviations, and 
correlations 18 

Multiple Linear Regression 

ORDER — rearrangement of intercorrelations 19 

MULTR — multiple regression and correlation 20 

Polynomial Regression 

GDATA — data generation 21 

Canonical Correlation 

CANOR — canonical correlation 22 

NROOT — Eigenvalues and Eigenvectors of a special, 
nonsymmetric matrix 23 



Analysis of Variance 

AVDAT — data storage allocation 24 

AVCAL — Sand A operation 25 

MEANQ — mean square operation . 26 

Discriminant Analysis 

DMATX — means and dispersion matrix 27 

DISCR — discriminant functions 28 

Factor Analysis 

TRACE — cumulative percentage of Eigenvalues 29 

LOAD — factor loading . 30 

VARMX — varimax rotation 31 

Time Series 

AUTO — auto co variances 32 

CROSS — crosscovariances . . . 33 

SMO — application of filter coefficients (weights) 34 

EXSMO — triple exponential smoothing 35 

Nonparametric Statistics 

2 

CHISQ — X test for a contingency table 36 

UTEST — Mann-Whitney U test . 37 

TWOAV — Friedman two-way analysis of variance 38 

QTEST — Cochran Q test 39 

SRANK — Spearman rank correlation 40 

KRANK — Kendall rank correlation 41 

WTEST — Kendall coefficient of concordance 42 

RANK — rank observations 43 

TIE — calculation of ties in ranked observations 44 

Random Number Generators 

RANDU — uniform random numbers 45 

GAUSS — normal random numbers 46 



Special Matrix Operations 

MINV — matrix inversion. 47 

EIGEN — Eigenvalues and Eigenvectors of a real, 

symmetric matrix 48 

SIMQ — solution of simultaneous linear, algebraic 

equations 49 



Matrices 



GMADD — add two general matrices 50 

GMSUB — subtract two general matrices 51 

GMPRD — product of two general matrices 52 

GMTRA — transpose a general matrix 53 

GTPRD — transpose product of two general matrices ..... 54 

MADD — add two matrices 55 

MSUB — subtract two matrices 56 

MPRD — matrix product (row into column) 57 

MTRA — transpose a matrix 58 

TPRD — transpose product 59 

MATA — transpose product of matrix by itself 60 

SADD — add scalar to matrix 61 

SSUB — subtract scalar from a matrix 62 

SMPY — matrix multiplied by a scalar 63 

SDIV — matrix divided by scalar 64 

RADD — add row of one matrix to row of another matrix .... 65 
CADD — add column of one matrix to column of another 

matrix 66 

SRMA — scalar multiply row and add to another row 67 

SCMA — scalar multiply column and add to another column. . . 68 

RINT — interchange two rows 69 

CINT — interchange two columns 70 

RSUM — sum the rows of a matrix 71 

CSUM — sum the columns of a matrix 72 

RTAB — tabulate the rows of a matrix 73 

CTAB — tabulate the columns of a matrix 74 

RSRT — sort matrix rows 75 

CSRT — sort matrix columns 76 

RCUT — partition row-wise 77 

CCUT — partition column-wise , 78 

RTIE — adjoin two matrices row-wise 79 

CTIE — adjoin two matrices column-wise 80 

MCPY — matrix copy 81 

XCPY — copy submatrix from given matrix 82 



RCPY — copy row of matrix into vector 83 

CCPY — copy column of matrix into vector 84 

DCPY — copy diagonal of matrix into vector 85 

SCLA — matrix clear and add scalar 86 

DC LA — replace diagonal with scalar 87 

MSTR — storage conversion 88 

MFUN — matrix transformation by a function 89 

RECP — reciprocal function for MFUN 90 

LOC — location in compressed-stored matrix 91 

ARRAY — vector storage-double dimensioned storage 

conversion , 92 



Other Mathematics 



Integration 



QUADR — integral of tabulated function by quadrature .... 93 
SMPSN — integral of given function by Simpson's rule .... 94 
RK1 — integral of first -order differential equation by 

Runge-Kutta method 95 

RK2 — tabulated integral of first -order differential 

equation by Runge-Kutta method 96 

RK3 — tabulated integral of a system of six first -order 

differential equations by Runge-Kutta method 97 

Fourier Analysis 

FORIF — Fourier analysis of a given function 98 

FORIT — Fourier analysis of a tabulated function 99 

Special Operations and Mathematical Functions 

GAMMA — gamma function 100 

LEGEN — Legendre polynomial 101 

BESJ — J Bessel function 102 

BESY - Y Bessel function. 103 

BESI — I Bessel function 104 

BESK - K Bessel function 105 

CELI1 — elliptic integral of the first kind 106 

CELI2 — elliptic integral of the second kind . . 107 

EXPI — exponential integral 108 



SICI — sine and cosine integral 109 

CS — Fresnel integrals HO 

Roots of Nonlinear Equations 

RTWIT — refine estimate of root by Wegstein's iteration ... Ill 
RTMIT — determine root within a range by Mueller's 

iteration . , . 112 

RTNIT — refine estimate of root by Newton's iteration .... 113 

Roots of Polynomial 

POLRT — real and complex roots of polynomial 114 

Polynomial Operations 

PADD — add two polynomials 115 

PADDM — multiply polynomial by constant and add to 

another polynomial 116 

PC LA — replace one polynomial by another 117 

PSUB — subtract one polynomial from another H8 

PMPY — multiply two polynomials . . 119 

PDIV — divide one polynomial by another , . 120 

PQSD — quadratic synthetic division of a polynomial 121 

PVAL — value of a polynomial 122 

PVSUB — substitute variable of polynomial by another 

polynomial » , 123 

PCLD — complete linear division 124 

PILD — evaluate polynomial and its first derivative ..... 125 

PDER — derivative of a polynomial 126 

PINT — integral of a polynomial 127 

PGCD — greatest common divisor of two polynomials. .... 128 

PNORM — normalize coefficient vector of polynomial 129 



SUBROUTINE TALLY 



**<.****■****:£*** 

FNT1-* 

******** ***-$*#,!t 



* + *************** 

* *Nf) ** * * 

* * * * CALCULATE * 

AIUI LESS +*+**w*+#*CO***>* MEANS FOR * 

THAN VHIN(I) * ** * *EACH VARIABLE* 

* * # ** * * 

* * * *************** 



*************** 



************* 



REPLACE 

vninTii by 

AHJ) 



*************** 



*************** 
* CALCULATE * 



*************** 



*| ***************** 



*************** 

* * 
INITIAL I/F 

*VMlN AND VMAX 

* Vfr.TIRS 

* * 

!=*#*# ********** 



*************** 



*************** 



£*-£*«********** 

* * 

* 
*^***-*** ******* 



*************** 

* * 

* REPLACE * 
I* VMAXU) BY * 

* A(IJ> * 

* * 
*************** 



*YE5 *** 

* * * 
SUMO *» B7* 

* * * 
* *** 

* * 



*************** 

* * 

* ACCUMULATE * 

* SQUARES OF * 

* ELEMENTS * 

* * 
*************** 



*************** 



.*.;<************* 



********* 



*************** 

Vf * 

* SET 1=1 

* * 
*************** 



*************** 



♦INCREASE I BY****) A7* * 

* 1 * * * * 

* * *** * 
*************** * 



******* *********:!<** 



*************** 
** * * 

* * RET ELEMENT * 
■,--A7***>* MIJ» IN * 

* * tms'ERVATIUN * 
** * MATRIX * 

*************** 



*YES *** 

* * * 
*} CO* 

* * * 
* *** 



************ *** 



******* ****** ** 



* ACCUMULATE ****) BO* 

* TOTALS * * * 

* * *** 
*************** 



♦INCREASE J BY****I A4* 

* 1 * * * 

* * *** 
*************** 



SUBROUTINE BOUND 



** ****** ******* 

ENTER 
*************** 











* * 












• 


* 






*• . 




• 






♦NO 




* * 


* 








* 




•SO*) 








J»NO 




•••••*••• 


* * 


* 








* 


* 


** 




* 


* 


* 
• • 

♦ ves 


* 


* 
* 
* 
* 



.* 



*************** 

* * * 

* CLEAR OUTPUT * * EXIT 

J VECTORS ( * * „«,,„*,**,„„ 

*************** 



* 

+ 

* 

******************* 



* SET J=l * +INCREASE J 8V* 

* * 2 l I 

* * * * 
*************** *************** 



* ( ******************************** 
* 

*YES *** 

* * * 
| J ) = *l BO* 

+ * * 
* *** 



* set r-i 

* * 

*************** 



*!>*« ft ********** 

** * * 

* * * GET ELEMENT * 

A=i***)* At IJ» IN * 
:■ * * OBSERVATION * 

** * MATRIX * 

*************** 



*YES 



***************)* ADO I TO ***** A9* 

* UN0EMU * * * 

* * *** 
*************** 



* * *************** 

* *YES * * *** 

AIIJ1 GREATER * *♦****.*****»**>* ADO 1 TO J***l A9* 
THAN 8HKI1 * * OVERtll * \,* 

* * * * *************** 



* ADD 1 TO » 
» BETWUI * 

* * 
*************** 



*ND 



*************** 



*** 



I,NV " *************** >*1NCREASE I BV****» A5* 

* * * l * *** 

t *************** 

* * *** 

* YES * * 

**********)* BO* ! 

* * 
*•* 



SUBROUTINE SUSST 



*************** ** * *NQ 

* * * * * 

ENTER * •BO*J I-NO ********* 

* * * * * * 

*************** ** * * * 



***#tt ********** V 

s * *************** 



* * *************** 

*************** 



********** ********* 



\* SET J=l * *INCREASE I fly****) 42* 

* * * ! * * * 

* * * * *** 
*************** *************** 



*************** 



*************** 

* * 

* GET ELEMENT * 

* IN * 

* OBSERVATION * 

* MATRIX * 
*************** 



*************** 

* * 

* GET * 

* RFLATIONAL * 

* CODE IN C * 

* MATRIX * 
*************** 



CONDITION OF ***************!* SET RIJ) ■ 1 

C SATISFIED * * 

+ * * * 

* * *************** 



*(******************•*******•***** 



*************** 



*************** )*INCREASE J BY****1 A3* 

* 1 * * * 

* * •** 
*************** 



*************** 

*B * *** 

******* ******* * * 

* CALC SI I) ***** BO* 

* WITH tlOOL * * * 

* SUBR * *** 
*************** 



SUBROUTINE MSNT 



**** *********** 



*************** 



* SET J-l TO 

♦begin Testing 

* from first 

* ROM OF A » 

*************** 



*»2*»*)» S(J>-1 * 

*< * ■* * 

** * * 

*************** 



* *NO 

MISSING OR * 
ZERO VALUE ********* 
FOR J-TH ROW * * 



* SIJ1-0 * 

* * 

* * 
*************** 



*(***************** 

* 

V 
* 
* 

♦ YES 
* 

********* 



* AOO 1 TO J ****! A2* * 

* * * * * 

* * *** * 
*************** * 



******************* 



*************** 



SUBROUTINE TASL 



* * *************** 

*************** ** • *YES ** * ACCUMULATE - 



TOTALS AND 
a*********CO***|* SQUARES OF 



.*****«»***•*** ** * _ _ * * ** *,iEa^j«JK!»J,J 



* TOTALS 

* (EQUATION 
************ 



, * * ***** DETERMINE * 

. INITIALIZE * JINCREASE J BY***») A3* * * INTERVAL FOR* 

*VH[N ANO VMAX * * 1 * * * * * * IIJI I 

^ * * * *•* * * * 

*************** *************** * *************** 



******************* 



• *NO 

* * LONER LIMIT * 

SET J * l .* * , SAME LIMl¥ P,,,ER * •"***". * FREO VECTOR * 



*************** 



*************** 



♦ YES *** * * * • *Yti 

* * * *REPLACE LOWER* • * . „„ * »*,*,.«.. 
•| 80* * LIMIT BY * * * J » NO ********* 

* * * * VMIN, UPPER * * * * J 

* *** *LIMIT BY VMAX* * * „ * J 

*************** * * * * 



*(****•******•***** 

* 

+ 

V 

*************** 

* * * % * * *"% 

* '"'"(IjflN 1 I I CLEAR OUTPUT • .INCREASE J 8Y****I B8* 
« OBSERVATION * * AREAS * * 1 * '„„* 

t MiTDIV t ♦ * * »*■» 

*************** *************** - *»«***********« 



******************* 



* * CALCULATE * * * * CALCULATE * 

A(IJ) LESS ********* .INTERVAL SHE* *C5***I* RELATIVE * 

THAN WIN • * . I EQUATION!.. *„* S f ISS^f 1*85 E | «t 

* * • *****•****••*** *************** 



CALCULATE * 

MEANS * 

(EQUATION 41* 



*************** 



*g **************** 4 



*************** *************** 



*************** 

*YE£ *** 



, , * CALCULATE « 

AtU) GREATER *l BO* * S(J> - *l C5* * n |J*NDARO J 

THAN VHAX *** * *** * DEVIATIONS * 

I HAN VN»« ***** * * *** * IE0UATION 51* 

* * * * *************** 



*************** 

i ****** GET ELEMENT * * „ „ 

* REPLACE VMAX****I BO* *B8**«I* All J] IN * * EXIT 

. BYA.U) J ..„. *„* .OBSERVATION* *,***„.*,.*..„ 

*****••*••***** *•*»•**»****»** 



* ADO 1 TO •***» CO* 

* COUNT * * * 

* • *** 
*************** 



SUBROUTINE TAB2 



*************** 
ENTER 

*************** 



•*********C0*** 



*************** 
* VARIABLE • 



* *NO *** 
LOWER LIMIT • * * 

SAME AS UPPER *) B4* 

LIMIT • * * 

* * *** 



♦INCREASE J BV***»I A** » 

* I + * * * 

* * #** * 
*************** * 



****##***♦***** 

* GET SECOND * 

* VARI-- " 



*************** 



******************* 



* INITIALIZE 

* VMIN, VMAX 



*************** 



*************** 

* * 
♦REPLACE LOHER* 

♦ LIMIT BY * 

• VMIN, UPPER * 



*************** 

* DETERMINE • 

* INTERVAL FOR* 

* SECOND * 

* VARIABLE * 

* ELEMENT • 
*************** 



*************** 
** * * 

* * * 
«A3***>* SET J = 1 

* * * 

*************** 



*NO **♦ 
* * * 
SECOND ♦) A3* 

VARIABLE * • » 
* *** 



*************** 
ADD 1 TO * 
*CORR6SP|jND|NG* 

JENCV * 
TRIX * 
*************** 



* * 


45 V 


* * 


******* ******** 


** * *VES *** 


** * * 


* * # * * * 


* * * CALCULATE * 
JB^I.INTgRVAL^EJ 


*A4*t S(JI - *> BO* 

* # * * * * 


** * * *** 


** * * 



♦***********♦♦* 



••••*••**••**** 

* ACCUMULATE * 

* TOTALS ANO * 

* SQUARES OF * 

* TOTALS IN * 

* ST ATI * 
*************** 



*************** 

* * 

* GET ELEMENT * 

* AIIJI IN * 

* OBSERVATION * 



*************** 



*** ************ 



* CLEAR OUTPUT 

* * 
*************** 



*************** 

* ACCUMULATE * 

* IQTALS AND * 

* SQUARES OF * 

* TOTALS IN * 

* STAT2 * 
*************** 



Al Ul LESS 
THAN VMIN 



********* 



******* ******** 

* * 

* SET J » I 

* 

* * 
•*****••**•**** 



♦YES **** 

* * * 
*l AO * 

* * 2* 
* **** 



20 V 


* * * 


V 


*************** 


* * * 


*************** 


* * 


* ** * *YES *** 


* • *** 


* * 


** * * * * * 


* * * • 


* REPLACE VMIN* 


**BT*1 SCJI -0 *l C6* 

** * * * * * 


♦INCREASE J BY****j B7* 


* BY A(U) * 


* * 


* ** * , * . *** 


* * *** 


*************** 


* * * 


*************** 


+ 


* * * 




* 


* * NO 





• I ***************** 



*n 



************ 



AIIJI GREATER *l BO* 

♦ THAN VHAX * * ♦ 
* * *** 



A ?8UNT T ° 



* REPLACE VMAX***M BO* 
« BY AIIJI * * * 

* • •** 
*************** 




SUBROUTING TAB2 



100 

*************** 
*** * * 

* * * CALCULATE * 
«AC **l* KELATIVE * 

• * * F«EOUENCIES * 
««* * (EOUATION 2)* 

*************** 



* CALCULATE * 

* MEANS • 

* (EQUATION 3>* 

* * 
*************** 



* STANDARD * 

* OEVIATIONS * 

* (EOUATION 41* 
*************** 



EXIT 

*************** 



SUBROUTINE SUBNX 



******** ******* 



enter . *.<», M^flun.^n* 1 : 

*************** ** * "*" * *** 



* * *COUNT NONZERO* 

* SET L«0 LL-0 * * COOES IN * 

* * * VECTOR S * 

* * * * 

*************** 



FDR EACH *************** 

COLUMN OF * 

OATA MATRIX * EXIT 



A* UU 

FOLLOWING 



*************** 



*************** 



*************** 



*A4***t* ADD 1 TO L * 

* * * * 

** * * 

*************** 



********* 



* ADD 1 TO 11 * 

* * 

* * 
*************** 



* MOVE DATA * 

* FROM A(L> TO* 

* D<LL> * 

* * 
*************** 



*{***************** 



*YES *** 



* ADD 1 TO I ***** A4* 

* * * * 

* * *** 

*************** 15 



SUBROUTINE NONCN 



*************** ** * *NO 

• * * * * 
ENTER * *BO*> OPTION - * 

* * * * * 
*************** ** * * 



• CALC FOURTH * 
•MOMENT EQ.I6I* 

* * * * 

*************** ******•*•**♦*•* 



*i ***************** 



V 

* * ** ft************** 

* C»LC NUMBER * * * * 

* OF FREQUENCY* *S2*I* EXIT 

* INTERVALS * * * * 

* EQ. Ill* ** *************** 
*************** 



*************** 



* CALC SUM OF • 

* ALL * 

* FREQUENCIES * 

* Ea.i?i * 
*************** 



* CALCULATE * 

* FIRST MOMENT************ 

* Eg. (31 * * 

* * * 
*************** *- 



IF OPTION * 

5, CALC 

SECOND. 

THIRD, ANO 

FOURTH MOMENT 



******************* 



OPTION = 2 ********* 



* CALC SECOND ****) »Z* * 
•MOMENT EO.I4I* * * « 

* * •** * 
*************** * 



******************* 



* CALC THIRD ****) B2* 
•MOMENT EO.I*)* • * 

* * *** 
*************** 



16 



*************** 



SUBROUTINE TTSTT 



I* XnI^oIgRe'e's' 6 ************ 



OF FREEDOM 
*************** ** * EQ. (5)-(7»_ 



*************** 



(ALC.i^yHUE*,. 



INITIALIZE « *B1**+1* AND DEGREES *****»*****> 

WORKING ARFA * * * * OF FREEDOM * J. 

- 'i****** * 



*************** ******** 



*************** 

■* * ** * * 

* * * * * C ALC T-VALUE* 

* CALCULATE « *B2***|* AND DEGREES * 
*MEANS EQ. (1)* * * * OF FREEDOM * 

* * ** *EQ. UOl-lUI* 

*************** *************** 



*| ***************** 
* 
* 
200 



V 

* *YES ** *************** 
OPTION * * * * 

GREATER THAN **********83*»* EXIT 

4 * * * * 

* * ** *************** 



*YES *** 
* * * 
OPTION » t, *» B2* 



*YES *** 



OPTION LESS 

THAN 1 

• * 



* CALCULATE * 
•VARIANCES EO.* 

* 121 * 

* * 
*************** 



********* 



* AND DEGREES ****! B3* 

* OF FREEDOM * * * 

* EO. (31-14) * *** 



*************** 



******************* 

* 
V 

* 170 

* 

•NO *** 

* * * 
ION » ? *) Bl* 

* * * 
* *** 

* 

* *** 



17 
*** 
160 



SUBROUTINE CORRF 



•****••*****«** 

ENTER 
*************** 



*************** 



*30***l* 5UH 
• ERt 



• OBSERVATIONS* 

* * 
*************** 



*************** 



I* SUBTRACT 1 * 

* FROM KK * 

* * 
*************** 



*************** 

* * 

* INITIALIZE 

* WORKING 

* STORAGE 

* * 
*************** 



*************** 



* MEANS FROM H* 
TIO. 



♦NO *•* 

* * * 
*) B5* 

* * * 
* **• 



ARE DATA 

ALREADY IN 

CORE 



********* 



*************** 

* * 
*SUM DEVIATION* 

* FROM * 

* TEMPORARY * 

* MEANS * 
*************** 



*************** 

* * 

* CALCULATE * 

* MEANS. SEE * 

* EQUATION (31* 

* * 
*************** 



*************** 

* * 

* CALCULATE * 

* MEANS. SEE * 

* EQUATION (31* 

* * 
*************** 



*************** 

* SUM CROSS * 

* PRODUCTS OF * 

* DEVIATIONS • 

* FROM TEMP. * 

* MEANS * 
*************** 



20S V 

*************** 

* CALC SUMS OF* 

* CROSS * 

* PRODUCTS OF * 

* DEVIATIONS. * 

* EQ (II * 
*************** 



*************** 

* * 

* SUM * 

« DEVIATIONS * 

* FROM MEANS * 

* * 
*************** 



*************** 



• CALCULATE * 

• KK-N-M * 

• * 
*************** 



*************** 

* CALC CORR. * 

* COEFF, IN * 

* SYMMETRIC * 

* MATRIX R. EO* 

* 1*1 * 
*************** 



*************** 

* CALC SUMS OF* 

* CROSS * 

* PRODUCTS OF * 

* DEV. FROM * 

* MEANS * 
*************** 



* *************** 
» ** *DATA * 
** * ******* ******* 

**B5*«*I* GET AN * 
** * * OBSERVATION * 

* *********••••** 



*************** 

* * 

* CALCULATE * 

* STANOARO * 

* DEVIATIONS, * 

* EQUATION (SI* 
*************** 



* 205 

*** 

PFRFORH * • 

ADJUSTMENT BY **) C3* 

EQUATION (II * * 

*** 



*************** 



* SUM THE * 

* OBSERVATION * 

* * 
*************** 



*************** 
•COPY DIAGONAL* 

* OF SUMS OF * 

* CROSS * 

* PRODUCTS * 

* MATRIX * 
*************** 



******************* 



127 V 

*************** 

♦DATA * 

******* ******* 

* GET AN * 

* OBSERVATION * 

* * 
*************** 



*•****••**•*•*• 

* SUBTRACT * 

* TEMPORARY • 

* MEANS FROM * 

* THE * 

* OBSERVATION * 
*************** 



*************** 



*************** 



•DATA" IS 

CALLEO M 

TIMES TO GET 

M 
OBSERVATIONS 



******•***•••** 

* * 
♦SUM DEVIATION* 

* FROM • 

* TEMPORARY • 

* MEANS * 
*************** 



*************** 

* * 

* CALC * 

* TEMPORARY ** 

* MEANS. SEE * 

* EQUATION (21* 
*************** 



*************** 

* SUM CROSS * **• 

* PRODUCTS OF * * • 

* DEVIATIONS *•••] CO* 

* FROM TEMP. « * • 

* MEANS * ••* 
*************** 



18 



SUBROUTINE ORDER 



*************** 



*************** 



*************** 

* CUPY CORR * 

* COEFF (X VS * 

* Yl IN RYIJI * 

* FOR RET. * 

* VARIABLES * 
*************** 



*************** 

* COPY * 

* INTER-CQRR * 

* COEFF OF * 
*HETAINFD INO.* 

* VARIABLES * 
*************** 



*************** 

* STORE * 

* DEPENDENT * 

* VARIABLE * 

* NUMBER IN * 

* ISAVE1K+1 I * 
*************** 



*************** 



*************** 



SUBROUTINE MULTR 



******** ******<! 

FNTE^ 
*************** 



*************** 
** * * 

* * » CALCULATE * 
♦BO***l*T-VALUES. SEE* 

* « "EQUATION (12)* 
** + * 

*************** 



*************** 



*************** 

* ■ * 

* CALCULATE * 
•STO. ERROR OF* 
•ESTIMATE. SEE* 
•EOUATION (101* 
*************** 



*************** 

* CALCULATE « 

* COEFF OF * 
*0 FTC'! MI, NATION* 

* . SEE * 

* FOUATION ( '*)* 
*************** 



*************** 

* * 

* CALCULATE * 

* F-VALUE. * 

* EOUATION I8»* 

* * 
*************** 



*************** 

* * 

* CALCULATE * 

* A EGRESSION * 

* COEFF. SFE * 

* FOUATION (21* 
*************** 



*************** 

•STORE RESULTS* 

* OF * 

* CALCULATION * 

* IN ANSU-IOI* 
*************** 



*************** 

* * 

* CALCULATE * 

* INTERCEPT. * 
« SEE EQUATION* 
<• HI * 
******** ******* 



*************** 

EXIT 
*************** 



*************** 

* * 

* CALCULATE * 

* SSAR. SEE * 

* EOUATMN (A)* 

w * 

*************** 



*************** 

* CALCULATE * 
■('MULTIPLE CORR* 

* CDEFF. SEE * 
« FOUATION (51* 

*************** 



*************** 

* + 

* CALCULATF * 

* SSDK. SEE * 

* EQUATION 171* 

* * 
*************** 



*************** 

* * 

* CALCULATE * 

* VARIANCE OF * 
*tSTIMATE. SEE* 

* FOUATION (9>* 
*************** 



*************** 

* CALCULATE * *** 

* STO. OEV. OF* * * 

* HFC. COFFF. ****) RO* 

* SEE EOUATION* * * 

* (11) * *** 

*************** 



SUBROUTINE GOAT A 



*************** 

ENTER 
************ **+ 



*************** 

* * 

* GENERATE X * 

* VARIABLES. * 

* EQUATION ID* 

* * 
*************** 



» CALCULATE * 

* MEANS. * 

* EQUATION (2)* 

*************** 



* CALC SUMS OF* 
•CROSS-PRODUCT* 
*OF DEVIATIONS* 

* EQ. (3) * 

*************** 



* CALCULATE * 

* CORRELATION * 

♦COEFFICIENTS.* 

* * 
*************** 



* DEVIATIONS. * 

♦ EO. <5) * 

*************** 



EXIT 

*************** 



SUBROUTINE CANQR 



*************** 



*************** 



NUMBER OF 
CRITERIA! 



********* 



*************** 

* * 

* PARTITION * 

* CORRELATION * 

* MATRIX * 

* EQUATION (1»* 
*************** 



*************** 



*************** 



******************* 



*************** 

* * 

* CALCULATE * 

* EIGENVALUES,* 

* EIGENVECTORS* 

* EQUATION (2)* 
*************** 



*************** 

* * *** 

* * * * 

* ADD 1 TO I ****) A5* 

* * * * 

* * *** 
*************** 



*************** 

*NROOT * 

******* ******* 

* EIGENVALUES * 

* EIGENVECTORS* 

* * 
*************** 



*************** 



SET I - 1 



*************** 



** 



*************** 



CALCULATE * 

*A5***)+ CANONICAL * 

* * * CORRELATION * 

** * EQUATION (51* 

*************** 



*************** 

* * 

* CALCULATE * 

* CHI-SQUARE * 

* EQUATION (6)* 

* * 
*************** 



************ *** 

* CALC DEGREES* 

* OF FREEDOM * 

* FOR * 

* CHI-SQUARE * 

* EQUATION (71* 
*************** 



*************** 

* * 

* CALCULATE * 

* RIGHT HANO * 

* COEFFICIENTS* 

* EQUATION I8»* 
*************** 



*************** 

* * *** 

* CALCULATE * * * 

* LEFT HAND ***** 80* 

* COEFFICIENTS* * * 

* EQUATION (91* *** 
*************** 



SUBROUTINE NROOT 



*************** 



*************** 



V 

*************** 
*FIGEN * 

******* ******* 

* EIGFNVALUFS,* 

* EIGENVECTOR * 

* * 
*************** 



V 

*************** 

* FORM RECIP. * 
*|]E SOUA«RDOT* 
*' OF * 

* FIGENVALUFS.* 

* EO. (1) * 
*************** 



* FORM * 

* a**(-i/2i * 

* * 
*************** 



* FORM 

* SYMMETRIC 
♦MATRIX " "- 



TRIX S. EO.» 
(2>-OI * 

************ 



V 

*************** 
*EIGEN * 

******* ******* 

* EIGENVALUES,* 
« EIGENVECTOR * 

* * 
*************** 



* NORMALIZE * 
•EIGENVECTORS.* 

* 60. UI-I5I * 

* * 
*************** 



EXIT 

*************** 



SUBROUTINE AVDAT s 



*************** 

ENTER 
** ************* 



*************** 

FXIT 

*************** 



*************** 
+ * 

* CALC SIZE OF* 
♦OATA ARRAY X-* 

* SEE EQUATION* 

* til * 
*************** 



*************** 



*************** 



*************** 

* * 

* MOVE DATA TO* 
♦UPPER PART OF* 

* ARRAY X * 

* * 
*************** 



KOUNT(J) * 

LEVEUJJ 



**:******* 



*************** 

* CALC STEPS * 
*USF_D TO STORF* 

* DATA. SEE * 

* EQUAT [ON * 

* 12-3) * 
*************** 



*************** 



*************** 



*************** 

* INITIALIZE * 

* KOUNT(J)=l * 

* FOR DATA * 

* SUBSCRIPTS * 

* FOR J»t,K * 
*************** 



********* 



*************** 



*************** 



* ADD 1 TO J ****! B2* ** 



*************** 



*************** 



******************* 



*************** 
*♦ * * 

* * *CALC LOCATION* 
*A6***»* (L) TO STORE* 

* * * A DATA. SEE * 
** * EQUATION (41* 

************ *** 



SET 

KQUNTIJI-l 



*************** 



*) B4* 
* * * 
*•* 



****************** 



*************** 

* * 

* REDISTRIBUTE* 

* I-TH DATA * 
♦POINT TO XIU* 

* * 
*************** 



*************** 



AOO 1 TO I ****) A6* 



***** ** ****** * * 



*************** 



*************** 



*NO 


*** 


* 


* * 


1=0 


*» fll* 


* 


* * 


* 


*** 


4 




* 


*** 


* YES 


* * 


********** 


I* BO* 




* * 




*** 



SUBROUTINE AVCAL 



*************** 

ENTER 

*************** 



*YES ** * *NO 

* * * * * 
**********C0*> FINISH ALL 

* ** * * FACTOR RUNS * 
* * ** * * 

* * » 

* * * 

* * YES 



********* 



*************** 

* * 

* CALCULATE * 

* CONTROL FOR * 
*END OF FACTOR* 

* RUN * 
*************** 



*************** 



* ADD 1 TO J ****) A8* 



*************** 



*************** 



*************** 



******************* 



******************* 



*************** 

* * 

* LASTSU>^L + 1* 

* WHFRE LIS* 

* LOCATION OF + 

* LAST DATA * 
******** **** *** 



*************** 



*************** 



*************** 



ADD L TO I ****) A6* 



*************** 



*************** 

* * 

* LASTSI I) = * 

* LASTSt 1-1) +* 
*ISTFP< I), FOR* 

* 1=2, K * 
*************** 



*************** 
** * * 

* * * CALCULATE * 
*B3***i* NN*X(LU-SUM* 

* * * AND STORE IT* 
** * IN X(LL) * 

*************** 



*************** 

** * CALC DATA * 

* * * RANGE * 

*C3***>* CONTROL. C ** 

* L i- LASTS}!)* 



*************** 

* * 

* ISTEP( I) IS * 

* I-TH STEP * 
♦USED TO STORE* 

* DATA * 
*************** 



*************** 



ADD ISTEPID* 

TO LL * 



*************** 



C GREATER 
THAN 



********* 



*************** 



*************** 



********* 



*************** 



♦REPLACE L AND****) A7* 

* LL BY C+l * * * 

* * *** 
*************** 



******************* 



*************** 
** * * 

* * * SET L»l f 
*A6***)* LL*l L SUM=>0, 

* * * NN=LEVFL( I ) 

** * * 

*************** 



*************** 



ADO 1 TO J ***** 63* 



*************** 



*************** 



* ADD ISTEPII >****} A7* 

* TO L ANO LL * * * 

* * *** 
*************** 



******************* 



*************** 



*************** 



* * 167 

* * 

*N0 *** 

* * * 
L'LASTSm *J C3* 

* * * 
* *** 

* * 

* * *** 
* YES * * 
**********}* co* 



*************** 



******** ******* 



*************** 



* ADD ISTEPCU****) BO* 

* TO L * * * 

* * *** 
*************** 



SUBROUTINE HEANG 



*************** 

,**. t CALC FIRST * 

jboj«. : ANO y ||CaN D ; 

*************** ** :.*„„*„„*,: 



PERFORM MEAN ni5f5no°U 

SQUARE T SJ V nlrSi:Jl 

OPFRATION AS THE DEGREES 

FOLLOWS 0F FREEDOM 



* INITIALIZE * * OIWIDE I-TH J 

* WORK AREA AND * * SUM BY FIRST* 

* CONTROL FOR * * DIVISDR * 

* LAST DATA * * „„,„,„„ 
*************** *************** 



*************** 



*************** 



THE ABOVE 

RESULT IS SUM 

OF SQUARES 



******** ******* 

**** IcALC POSITION? LfiL^IS'UfcSKn*' 

„«»**>* IN OUTPUT * ! SUM n?X.i£E° ND I 

t * * AREA TO SUM * * DIVISOR * 

** * NN-TH DATA * * „„„„„,! 

*************** *************** 



* SQUARE NN-TH* IB!,,??°YI 

* DATA AND * RE SULT IS 
♦ACCUMULATE IN* MEAN SQUARES 
*THE ABOVE LOC* 

*************** 



«y ES ^ * .Nil 

. LAST * ********* * I - N ********* 

DATA ♦ * * « * 



*************** * 

, * *«* * 

* * * * * 

* ADO 1 TO NN ****) A4» * * EXIT 

J * *** * *************** 

*************** * 



******************* 



NUMBER OF * J J * * 

OUTPUT N Is 2 * * ADD I TO I ****» BO* 

TO THF K-TH 2 1 * *** 

™WCK * *************** 



******************* 



*************** 



*************** 



SUBROUTINE DHATX 



*************** 
****** ********* 



* * * CALC POOLED * 
*BO***)* DISPERSION • 



MATRIX. 

*************** *• * EQUATION (31* 

* *************** 



* INITIALIZE * * __„. 

* WORKING * * EXIT 

* STORAGE t * * *„„„„.„.„ 
*************** 



*************** 



*************** 

** * SUN K-TH * 

* * * GROUP DATA * 

4A-i***(*AND CALCULATE* 

•i 1 * MEANS * 

** * EQUATION <1»* 

*************** 



* SUBTRACT * 

* MEANS FROM * 

* DATA * 

* * 
*************** 



* CALC SUM OF * 
*CROSS PRODUCT* 
*IJF DEVIATIONS* 

* EQUATION I 21* 
*************** 



*************** 

* * 

* ACCUMULATE * 

* SUM OF CROSS* 

* PRODUCT OF * 

* DEVIATIONS * 
*************** 



K - LAST 
GROUP 



*************** 

* * *** 

* * * * 

* ADD 1 TO K ****) A3* 

* * * * 

* * •** 
*************** 



ay 



SUBROUTINE DISCR 



*************** 



*************** 



*************** 

«* *FOR EACH CASE* 

* * *IN EACH GROUP* 

*BO***>* DO THE * 

« * « FOLLOWING * 

«* * STEPS « 

*************** 



*************** 

* * 

* CALCULATE * 

* COMMON MEANS* 
♦FOR VARIABLES* 

* EQUATION ( D* 
*************** 



*************** 

* * 

* CALCULATE * 

* DISCRIMINANT* 

* FUNCTION * 

* EOJATION 151* 
*************** 



*************** 

* CALCULATE * 

* GENERAL * 
♦MAHALANUBIS D* 

* SOUARE. * 

* EQ.I2) * 
*************** 



*************** 

* * 

* FIND THE * 

* LARGEST * 

* DISCRIMINANT* 

* FUNCTION * 
*************** 



*************** 



*************** 



*************** 

* CALCULATE * 
*PROBALITV FOR* 

* LARGEST * 
« FUNCTION * 

* EQ. til * 
*************** 



*************** 



SET K » 1 



*************** 



*************** 

EXIT 
*************** 



*************** 
** * CALC * 

* * * COEFFICIENTS* 
*A5***I* OF DISCRIMI-* 

* * *NANT FUNCTION* 
«• * EQ.O) * 

*************** 



*************** 



* CALCULATE * 

* CONSTANT * 

* EQUATION I'll* 

* * 

*************** 



•YES *** 



K » LAST 

FUNCTION 



*************** 



AOO I TO K ****) A5* 



*************** 



SUBROUTINE TRACE 



*************** 



*************** 



*************** 



• NO 
I-TH * 
EIGENVALUE ********* 



I ESS THAN 
CONSTANT 



*************** 

* * 

* CALCULATE * 

* CUMULATIVE * 
•PERCENTAGE OF* 
*K EIGENVALUES* 
*************** 



EXIT 
*************** 



******************* 



* ADD 1 TO K * 

* * 

* * 
*************** 



* DIVIDE I-TH * 
♦EIGENVALUE BY* 

* M. * 

* * 
*************** 



*************** 

* * *** 

* * * * 
« ADD 1 TO I ***♦! A2* 

* * * * 

* * *•* 
*************** 



SUBROUTINE LOAD 



*************** 

ENTER 
*************** 



************ *** 



*************** 



*************** 
** * CALCULATE * 

* * *SQUAREROOT OF* 
*A2***»* I-TH * 

* * * EIGENVALUE, * 
** * CALL IT SQ * 

*************** 



♦MULTIPLY I-TH* 

* EIGENVECTOR * 

* BY SQ * 

* * 
*************** 



********* 



*************** 



******************* 



* ADD I TO I ****» A2* 

* * * * 

* * *** 
*************** 
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SUBROUTINE VARMX 



*************** 

ENTER 
*************** 



** 



*************** 
* * 



** 



*************** 



*BO***l* ADO 1 TO NC * 

* * * * 

** * * 

*************** 



* » PERFORM . * 
*CO***l*OPERATI0NS IN* 
♦SECTION 18-31* 
* OF EQUATION * 
*************** 



** 



*************** 
•* * CALC * 

* * *OIFFER6NCE IN* 
* 0**«|*COMHUNALITIfcS* 

* * * SEE EQUATION* 
** * 1301 • 

**••******•*•*• 



v 

*************** 

* * 

* SET EPX - 
♦0.00116, LL » 

* K-l. NV » 1 

* * 
*************** 



NC GREATER 
THAN 3 



430 

*YES *** 
* * * 
*l C8* 

* * 
*** 



V 

*************** 
♦CALCULATE COS* 
* THETA. SIN * 

THETA. SEE * 



** 



*C 1***1* 
* * * 

** * iwt-uil ' 

*************** 



EQS. 
C141-U7!- 



*************** 
EXIT 

*************** 



*************** 

* * 

* SET NC * 0, 

* CONS - 

* 0.7071046 

* * 
*************** 



*************** 



** 

* * * 
*B2***I* SET J 



** 



*************** 

* 



*************** 

* * 
♦CALCULATE COS* 
*PHI, SIN PHI.* 
•SEE EQUATIONS* 

* C18I-I25) * 
*************** 



*************** 

* * 

* CALCULATE * 
SCOMMUNALITIES* 

* SEE EQUATION* 

* (21 * 
*************** 



** 



*************** 



*B3***)* SET Kl - J*l 
* * * 
** * * 

*************** 

* 



*************** 

* * 

* PERFORM * 
♦ROTATION. SEE* 

* EQUATIONS * 
♦126) ANO (271* 
*************** 



*************** 

* CALCULATE * 

* NORMALIZED * 

* FACTOR * 

* MATRIX. SEE * 

* EQUATION (31* 
*************** 



V 
*************** 

**** *CALCULATE NUM* 

*B4***t* AND DEN. SEE* 

* * * EQUATION <7>* 

** * * 

*************** 



*************** 



*C4***)* ADD 1 TO Kl * 

* . * * * 

** * * 

*************** 



*************** 
** * * 

• * * * 
eA 1 !***!* ADO 1 TO NV * 

* * * * 
** * * 

*************** 



*YES 



NUN LESS THAN 
DEN 

* * 



********* 



* I B4* 



*************** 

* CALCULATE * 

* VARIANCE FOR* 

* FACTOR * 

* MATRIX. SEE * 

* EQ. (4) * 
*************** 



•YES 



NUM GREATER 
THAN 06N 



********* 



*************** 

* * 

* * 

* ADD 1 TO J * 

* * 

* ♦ 
*************** 



NV GREATER OR 
EQUAL TO 50 



430 

•YES *** 
* * * 
*> C8* 
* * 
*** 



*N0 



V 
***•****♦♦***** 

* * 

* PERFORM * 

* CONVERGENCE * 

* TEST. SEE * 

* EQUATION (5)* 
*************** 



NUM ♦ DEN 
LESS THAN EPX 



250** 

** 

*** ** 

* * ** 

*> CO* ** 



* YES * * ** 
*********•)* C4* ** 

■ * * ** 

*** * 

******************* 

* 



* * *** 
* YES * * 

**********)* A5* 



430 

*************** 



*************** 

*, dhcoh * ***** * *"* * DENORMALliE 

*0PeSI?i8n§ InJ**.T Cl» * *C8***)« ..._FAC?oA 



T UrLl*M( iuiij -a ■ « - - 

♦SECTION (B-ll* 

* OF EQUATION * 
*************** 



~ MATRIX. SEE * 
•EQUATION (28)* 

*************** 



****************** 



*NO 



TEST 
SATISFIED 



*) BZ* 



* *•* 

* YES * * 

****•*****)* BO* 



320 V 

*************** 

* * *** 

.operators .n:..*t «: 

*************** 



*************** 

* CALCULATE « *** 

* FINAL * * * 
•COMMUNALlTIES****! DO* 
t ^E RATION* ♦„.* 

*************** 



SUBROUTINE AUTO 



*************** 

FNTFR 
*************** 



* 

* CAICUI.A II 

*riHF SPRITS A 

* 

t$**fr********* 



*#.***#******#* 



a************* f 



****:**#******** 

* l.r.LCM -Ml' * 
4 k'll.l.Y'VMU- * 



fr***iji*ii=+*4 
■XI r 
4********4 



32 



SUBROUTINE CROSS 



*************** 
ENTEtt 

*************** 



* GALCULATl- * 

* AVIFRAGfS HP * 

*T (ML SFSTFS A* 

* AN") ft * 

*************** 



** ************* 

* * 

* MULTIPLY THE* 

* :Wi ;i <AG£ (IF A* 

* HY THi- * 

* iW'IUGF flF H* 
*************** 



t'.rtl CtJLA'P." 

;,i»il!;sVA!'.I" 

AKCFS HP !.i 

LAGS A 

************ 



*************** 

* * 

* CALCULATE * 

* CMtSSVAfU- * 

* 5>m:es MP ?■; * 

* LiAiJS A * 

**S-**L ********** 



i***j}i******* *** 



*************** 
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SUBROUTINE SHO 



*************** 
ENTER 

*************** 



* CLEAR DUTPUT* 

* VECTOR R * 



*************** 



V 

*************** 

* CALC LOWER * 

* AND UPPFR * 
♦POSITIONS FOR* 

* SMOOTHED * 

* SERIES R * 
*************** 



♦SMOOTH SFR1FS* 
* A 8V WEIGHTS* 



*************** 



*************** 



*************** 



SUBROUTING EXSHO 



******** ******* 

FNTFR 
*************** 



* *YES 

* ARC * 
COEFFICIENTS ********* 

* A, B, C GIVEN * * 



* CALCULATE Ai* * 

* 9 c. ***********) 

* EQ.(5)-(7I * * 

* * * 
*************** * 



OH THE 
FULLflwINC, FOR 
At.l DATA* I = 



******************* 



*************** 

* * 

* FnRFCAST "NF* 

* PEKinO AHEAD* 

* FO.ll) * 

* * 
*************** 



*************** 

* * 
*UPD»TF OEFF.t 

hisinh current* 

* A, B, C * 
<• na. I 2 1-14 I * 

*************** 



*************** 

EXIT 
*************** 



36 



SUBROUTINE ClllSg 



*************** + * * 

* * * * 

ENTER * *BO»**)* SET ERROR 

* * « * CODE » 2 
*************** 



*************** 



V BO V 

*************** *************** 

* INITIALIZE * ** * COHPUTE * 

* ERROR COOE * * * * TOTALS OF * 

* AND * *B1***1* COLUMNS AND * 
*CHI-SQUARE TO * * * * ROUS IEO. * 

* ZERO * ** * 2-3) * 
*************** *************** 



*************** 

* * * * 

* COMPUTE * *COMPUTE GRANO* 
■* DEGREES OF * * TOTAL * 

* FREEDOM. * * (EQUATION *)* 

* (EQUATION II* * * 
*************** *************** 



********* * 2 BY 2 ********* 

« « SPECIAL CASE * * 



* 140 V 

* *************** 

* * * « COMPUTE * 

* « * * CHI-SOUARE * 

* SET ERROR * * * FUR OTHER * 

* CODE » } * * * CASES. * 

* * * * (EOUATION 6)* 
*************** * *************** 



*************** 



*************** * *************** 



******************* ******************* 

* * 

V * 

* 120 V 

* *************** 

*N0 * COMPUTE » 

* CHI-SOUARE * 

* FOR 2 BY 2 * 

* TABLE. • 

* (EOUATION 51* 
*************** 



*************** * V 

* * * *************** 



* * * *************** 

*************** * 



*************** 



******************* 



* *N0 

1/5 OF CELLS * 
CONTAIN LESS 
THAN FIVE * 



*** 

50 



SUBROUTINE UTEST 



*************** 



*************** 



* 

*ao* 



*************** 

*T I E * 

******* ******* 
**l* CONFUTE SUM » 

* OF T. teo.3>* 

* . • 
*************** 



*************** 
♦RANK * 

******* ******* 

* RANK SCORES * 

* fHOH TWO GRP* 

* * 
*************** 



********* 



SUM RANKS IN* 
LARGER GROUP* 



*************** 



Bgy! *I Q r;,* 



RD ************ 



******************* 



*************** 

* * 

* CALCULATE U * 

* PRIME, * 

* (EQUATION II* 

* * 
**** ****** ** ** * 



* COMPUTE 

* STANDARD 

* DEVI"'""' 

* (EQUA 
******** 



I MnUMKU -* 

VI ATI ON. * 
UATION 51* 
********** 



*************** 



*( **************** 



CALCULATE U.* 
(EQUATION ?)* 



*************** 



CALCULATE I.* 
(EQUATION 6>* 



*************** 



*************** 



*************** 



SET U - U 
PRIMR 



*************** 



*(***************** 



* *fJQ *** 

* * * 

N2 LESS THAN *) 80* 



*************** 



SET Z = 



*************** 



BO * 

V 
** ************* 

EXIT 

*************** 



SUBROUTINE THOAV 



*************** 

*************** ** * * 

* * * CALCULATE SUM* 
ENTER * *BO***>* OF RANKS * 

• * * * (EQUATION 11* 
*************** ** * * 

* *************** 



* * *************** 

* *YES * * 

* * ACCUMULATE * 
DATA RANKED ********* * SUM OF * 

* * * SQUARES OF » 

* * « *SUMS OF RANKS* 

* * * •*********••**• 



* * * LAST COLUMN 

* * * 

* * * 

*************** * * * 



* * * * CALCULATE * 
I* GET ROM IN * * * FRIEDMAN * 

* MATRIX A * * * STATISTIC * 
^ *«.*!* a ^ ^ t IE0UATI0N2I * 
*************** * *«*«*♦.******** 



* * CALCULATE * 

* * DEGREES OF * 

* * FREEDOM • 

* * (E0UATI0N3I * 

*************** * *************** 



*************** * 

♦RANK * * 

******* ******* * * 

* RANK DATA IN* * * EXIT 

* ROW * * * 

« * * *************** 

*************** * 



* MOVE RANKED * 
*ROW TO MATRIX* 



*************** 



*(*******•********* 



*************** 



* INITIALIZE 

* COLUMN INDEX 



*************** 



*A9***I*GET COLUMN IN****) BO* 

* * * MATRIX R * * * 

** * * *•* 

*************** 



SUBROUTINE OTEST 



*************** 



*************** 



*************** 



*BO***)* GET A COLUMN* 



*************** 



*************** 



* CLEAR TOTALS 



*************** 



*************** 

* * 
♦COMPUTE TOTAL* 

* OF COLUMN * 

* I EQUATIONS) * 

* * 
*************** 



*************** 



SET t » 1 



*************** 



* SQUARES OH 
♦COLUMN TOTALS* 

* * 
*************** 



*************** 



*A3***)* RET A ROW * 



********* 



*************** 



*************** 

* * 
♦COMPUTE TOTAL* 

* OF ROW * 

* (EQUATION 1>* 

* * 
*************** 



*************** 



♦INCREASE J 8Y****> BO* * 

* 1 * * * * 

* * *** * 
*************** * 



******************* 



* ACCUMULATE * 

* SQUARES OF * 

* ROW TOTALS * 

* * 
*************** 



*************** 



STATISTIC * 
■— ATIDN 3>* 



*************** 



** ** **** ******* 



*************** 

* * 

* FIND DEGREES* 

* OF FREEDOM * 

* (EQUATION *>* 

* * 
*************** 



********* 



*************** 



*************** 



*************** 



•INCREASE I BY****» A3* 

* l * * * 

* * *** 
*************** 



******************* 



*************** 



SET J » I 



*************** 



*| SO* 
• * * 

**♦ 



SUBROUTINE SRANK 



*************** ** * +N0 

* * * * * 

FNTFR * *BO*J TSA - ********* 

* * * * * * 

*************** ** * * * 



* *YES * *NO * 

* * * * 

DMA RANKED ********* * TS8 « ********) 



*************** * 

♦RANK * * 

******* ******* * 

* RANK DATA IN* * 

* VECTOR A * * 

* * * * * 
*************** * *************** 



******************* 



*************** * 

♦RANK * * * * 

******* ******* * * COMPUTE RS * 

* RANK DATA IN************ * C EQUATIONS * 

* VECTOR R ♦ ** * 4,5,61 *. 

* * ** * * 
*************** ** *************** 



******************* *(**************** 



*************** * * * 

■* * * * *NO 

* HOVE RANKED * * * * 

* DATA TG * * * N * 10 OR ********* 

* VECTOR R * * * LARGER * * 

* * * * * * 
*************** * * * * 



*( **************** 



*************** 

* * * * 

* INITIALIZE * * * 

* INDEX AT * * COMPUTE T * 
♦FIRST ELEMENT * * (EQUATION 7J* 

9 * * * 

*************** *************** 



*< ***************** 



*************** 

* * * * 

* FIND DEGREES* 

* OF FREEDOM * 

* (EQUATION 81* 

* * * * 
*************** *************** 



50 V 

*************** 

* ACCUMULATE * 

* SUM OF * * 

* SQUARES OF ♦ * EXIT 

* DIFFERENCES ♦ * 

* ( EQ. 1 > * *************** 
*************** 



*♦♦ 

* * 
1 A6* 

• * 
*** 



* * *** 

* COMPUTE TSA,^ ♦ * 
*TSB (EQUATION****! 80* 

♦ 2 1 * * * 

♦ ♦ *** 
*************** 



40 



SUBROUTINE KMNK 



*************** 



*************** 



*************** 



»BO***l* SST L « I, 

* * * » N*L 



*************** 



*************** 
** *TIE * 

* * ******* ******* 

*CO***>* COMPUTE TA * 

* * « (EQUATION II* 
** * * 

*************** 



* *YES ** 

* * * 

DATA RANKED **********B1* I 



G S!1J E 



********* 



*************** 

*TIE * 

******* ******* 

* COMPUTE TB * 

* (EQUATION D* 

* * 
*************** 



*************** 

*RANK * 

******* ******* 

* RANK DATA IN* 

* VECTOR A * 

* * 
*************** 



*************** 



* ADD I TO S * 



*************** 



TA AND TB - 



********* 



*( ***************** 



*************** * 

*KANK * * 

******* ******* * 

* RANK DATA IN************ 

* VECTOR B * 



*************** 



** 
** 



R(K) LESS 
THAN RUI 



********* 



* COMPUTE TAU ************ 

* (EQUATION 3)* ** 



*************** 



******************* 



******************* 



* MOVE RANKED * 

* DATA TO * 

* VECTOR R * 

* * 
**************# 



***** ********** 



SUBTRACT 1 

FROM S 



*************** 



* COMPUTE TAU * 

* (EQUATION 2»* 



*************** 



*(**************** 



*( ***************** 



* ( **************** 



♦SORT VECTOR R* 



*************** 



********* 



N LESS THAN 
10 



********* 



*************** 



* COMPUTE S * 



*************** 



******** ******* 



* ADD 1 TO L ***** Bl* * 



*************** 



90 * 
*************** 

EXIT 
*************** 



******************* 



******************* 



*************** 



*************** 



* COMPUTE SD, * 

* STANDARD * 

* DEVIATION * 

* (EQUATION 4)* 
*************** 



*************** 



*************** 



*************** 



* ADO I TO I * 



*************** 



*************** 



*************** 



*(******************************** 



*************** 



* GET RANK OF ****» 60* 

* VARIABLE B * * * 

* * *** 
*************** 



*************** 

EXIT 
*************** 



SUBROUTINE WTEST 



*************** 

ENTER 
*************** 



*************** 



*BO*+*l* SET J 



*************** 



* FREEDOM ♦ 

* (EQUATION 7J* 
*************** 



********#****** 



*************** 



*************** 



*B1***|* GET A COLUMN* 



*************** 



*************** 

EXIT 
*************** 



*************** 



GET A ROW 
VARIABLE 



*************** 



*************** 

* * 
♦CALCULATE SUM* 

* OF RANKS * 

* (EQUATION 2)* 

* * 
*************** 



* *NO 

* 
DATA RANKED ********* 



********* 



*************** 



*************** 



* ADO I TO J ****> Bl* * 



*************** 



*************** 



******************* 



******************* 



*************** 

♦RANK * 

******* ******* 

* RANK. DATA IN* 

* ROW * 

* * 
*************** 



* COMPUTE MEAN* 

* OF SUMS * 

* (EQUATION 31* 

* * 
*************** 



*{**************** 



*************** 

* TIE * 
******* ******* 

* COMPUTE T * 

* (EQUATION 1»* 

* * 

*************** 



♦ COMPUTE SUM * 
*OF SQUARES OF* 



******** **** *** 

* * 

* MOVE RANKED * 

* DATA TO * 

* MATRIX R * 

* * 
*************** 



*************** 



* COMPUTE W * 

* (EQUATION 51* 



*************** 



*YES *** 

* * * 
*l BO* 

* * * 
* *** 



*************** 

* * *** 

* * * * 

* ADO 1 TO I ****! AZ* 

* * * * 

* * *** 
*************** 



65 V 70 
*************** 

* * *** 

* COMPUTE * * * 

* CHI-SQUARE ****> CO* 

* (EQUATION 6)* * * 

* * *** 
*************** 



SUBROUTINE RANK 



*************** 

ENTER 

*************** 



*************** 

* * 

* INITIALIZE 
♦CONSTANTS AND 

* INDICATOR 

* * 
*************** 



*************** 

* * 

* FIND EACH * 
I* SUCCESSIVE * 

* SMALLEST * 

* ELEMENT * 
*************** 



* COUNT NUMBER* 

* OF SMALLER * 

* DATA POINTS * 

* * 
*************** 



* COUNT NUMBER* 
*OF EOUAL DATA* 

* POINTS * 

* * 
*************** 



IS THERE A ********* 
TIE * * 



* STORE RANK ************ 

* * ** 

* * ** 
*************** ** 



******************* 



* CALCULATE * 

* AVERAGE RANK* 

* AND STORE * 

* * 
*************** 



* < **************** 



LAST ELEMENT *l A2* 



EXIT 
*************** 



SUBROUTINE TIE 



*************** 

*************** *♦ * ACCUMULATE * 

* * * * CORRECTION * 
FNTFR * *BO***l+ FACTOR * 

* * * *tCT*CT*CT-CT»* 

*************** ** * /12- * 

* *************** 



* INITIAL IZF * * + * 
♦CONSTANTS AND * #91***1* RfSET 

* INDICATOR * * * * INDICATOR 

* * ** * * 
*************** *************** 



*< ********************* *********** 



* FIND EACH * 

* SUCCESSIVF * 
♦SMALLEST RANK* 

* * 

*************** 



********* 



90 * 

******** ******* 



*************** 



******************* 



* LOCATE TIES * 

* * 

* * 
*************** 



*************** 



*************** 



*ND *** 

* * * 
ANY TIES *) 61* 

* * * 
* **• 

* * 

* * 
* YES 



* * 80 
# * 

*YES *** 

* * * 

USE FIRST *) BO* 

EQUATION * * * 

* *** 



75 V 

*************** 

* * •** 

* ACCUMULATE * * * 

* CORRECTION ****> Bl* 
+ FACTOR * * * 

* CT*<CT-l>/?-* *** 
*************** 



SUBROUTINE RANDU 



*************** 
********-«*** **A 



*************** 

* * 

* COMPUTF IY * 

* USING INPUT * 

* IX * 

* « 
*************** 



*************** 

* * 

* CONVERT * 

* RESULT TO * 

* FLOATING PT * 

* IN RANGE- 0,1* 
*************** 



*************** 



*************** 



45 



SUBROUTINE GAUSS 



*************** 

FNT^R 
*************** 



*************** 

* * 

* COMPUTF 12 * 

* UNIFORM * 

* HANOnM * 

* NUMBERS * 
*************** 



*************** 

* * 
+AOO \2 RANDOM* 

* NUHBFRS * 
*■ TOGETHER * 

* * 
*************** 



** ****** **** *** 



* lOMPUTE y * 

* (ryuATUN z>* 



*************** 



*************** 

EXIT 

******** ****** -5 



SUBROUTINE MINV 



*************** 



*************** 



K GREATER 
THAN LOO 



********.**•**** 

* * 

* * * MULTIPLY * 

"EftMINANT — 



**********CO*** J*DETEft 



- D*BIGA * 

* * 

*************** 



*************** 

* * 

* INITIALIZE 

* DETERMINANT 

* AT I 

* * 
*************** 



*************** 

* * 

* INTERCHANGE * 
♦ROWS LCKI AND* 



*************** 



*************** 

* * 
♦REPLACE PIVOT* 

* BV ITS ************ 

* RECIPROCAL * * 

* * * 
*************** * 



*( ***************** 



*************** 



*A2***(* SET K 



K GREATER 
THAN MOO 



********* 



*************** 



******************* 



*************** 

* * 

* INITIALIZE 

* L(KI» M(K) * 

* RIGA 

* * 
*************** 



* INTERCHANGE * 

* COLUMNS MOO* 

* AND K * 

* * 
*************** 



********* 



*l ***************** 



*************** 

* * 

* INITIALIZE 
♦ARRAY POINTER 

* AT AtKK) 

* * 
*************** 



*************** 



********* 



•INCREASE K 
* 1 



BY**** I A 2* 



*************** 



******************* 



**** **** ******* 



*Ab***l* GET ELEMENT 
* * * A(IJ) 



*************** 



* SET * 

* DETERMINANT * 

*EOUAL TO ZERO* 

* * 
*************** 



* INTERCHANGE * 
♦ROWS AND COLS* 

* ACCOROING TO* 
*L AND K VECTS* 

*************** 



AUJ) GREATER 
THAN BIGA 



********* 



*************** 



*************** 



*************** 



*************** 



******************* 



*************** 

* REPLACE BIGA* 

* WITH AUJ). * 

* L(K» WITH I,* 

* MOO WITH J * 
*************** 



* DIVIDE COL K* 

* BY -BIGA ♦ 

* (EXCEPT * 

* PIVOT) * 
*************** 



* { ***************** 

* 

V 

* * 

* * 

* *YES **♦ 

* * * 

LAST ELEMENT ♦» BO+ 

* * * 

* * *** 

* * 

• * 
* NO 



*************** 

* reduce mat * 

* aTiji - * 

*AUKt*A(KJ) +* 

* AUJ). EX. I* 

* OR J«K • 
*************** 



*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ****> AS* 

* ELEMENT * * * 

* * *•* 
*************** 



*************** 

* * *** 

* OIVIDE ROW K* * * 

* BY BIGA •***» CO* 

* (EXCEPT * * * 

* PIVOT) * *** 
*************** 



SUBROUTINE E1GEN 



*************** *************** 

*************** ** * * ** * * 



ENTER * *BO***>*SET INDICATOR * *CO***>*INCREA$E I BY****) B3* *O0*> MV * I ********* 

* * * * TO 1 * *** i ****** * * 

*************** ** * * ** * * *** ** + * * 

* *************** *************** * * * 



*************** 

S * * **** * * 

* COMPUTE SIN * * * * SOLVE * * _,„ * 

********* * AND COS. * *ci***|* EQUATIONS * * „ , ^HPJ-rr.n * 

* * EQUATIONS * * * * (13>-(15) * * EIGENVECTORS* 

* * (3 1-J7I * ** * * * * 

* *************** *************** *************** 



*(******# ********** 



* GENERATE * ** **** * « 

* IDENTITY * * * SET I - 1 * *C2* I N = N ********* * EXIT 

* MATRIX* * * * *** * * * 

* * * * * *** * * *************** 

*************** * *************** * * * 



*(***************** * 

* 

V 

* * 

* * 

** * *YES * * *** 

* COMPUTE * *** * * *** 
♦INITIAL NORM.* *B3*> I - L ********* MNCREASE H BY****> A9* 

* FGUATION (11* *** * * * I *** 

* * *** * ** * *** 
*************** * * * *************** 



******************* 



* *YES * * *YES 

* COMPUTE FINAL* * * * * * 

* NORM. * * I « M ******#*) * h » N-l ********* 

* EQUATION (21* * * * * * * 

* *- * * ** * * 
*************** - ■ ■ - - * --■* +- * * * 



* SOLVE * + * * * * 

* EQUATIONS * * INCREASE L BY****1 A8* 

* <8I,(9I * * * i * * * 
* * * * * * * *** 
*************** *************** * *************** 



* ******************* 

* * 

* V 

* 150 * * 

* * * 

** * * * * * * *mq 

* * * COMPUTE * *REPLACE A«ILI* * * * 
*A6***|* THRESHOLD - + * - Bt AIIH1 ** * * INOICATOR = 1 ********* 

* * * NORM /* *C* ** ** 
** * NO. ROWS * * ***** 

*************** *************** * * * * 



* *( ***************** 

* * 

* V 
V 115 * * 125 

*************** * * 

** * * * *YES *** * * *** 

* * * * * *** * *** 

*A7***I* SET L = I * * MV " I *l «* * w CLEAR *} A7* 

«** * * *** * INDICATOR * * * 

** * * * * *** * * *** 

*************** * * *************** 



* * ******************* 

* * * 

* * V 

V 120 V 160 * * 45 

*************** *************** * * 

** * * * * * *yES *** 

* * * * * SOLVE * * THRESHOLD * * * 
*A8***1* SET M - L+l * * EQUATIONS * * GREATER THAN *l A6* 

* * * * * U0>-(12> * * FINAL NORM * * * 
** * * * * * * *** 

*************** *************** * * 



*YES ++* * * *** 

* * * * * * * 
« N *> CI* * SORT ****» DO* 

♦ * * * EIGENVALUES * * * 
* *** * * *** 

* *************** 
* *** 
* NO * * 
*•****••**)* CO* 





*YES 


*** 


AILNJ LESS * 


* * 


THAN 




*! C2* 


THRESHOLD 


* 


* * 




* 


*** 


* * 






* * 




**• 


* NO 




* * 


**********}* bo* 



SUBROUTINE SIMQ 



*************** 

ENTER 
*************** 



*************** 

** * * 

* * * eliminate * 
•b0***>* element from* 

* * * III other • 

** * ROUS * 

*************** 



*************** 

* * 
* 

* SET KS TO 

* 

* * 
*************** 



LAST COLUMN 



*************** 



* INITIALIZE 

* COLUMN tNDEX 



*************** 



*62***)*GET PRECEDING* 
* * * ROW * 



******* ******** 



*************** 



*A3***>* GET A COLUMN* 



*************** 



*************** 

* SOLVE EOU BY* 

* SUBSTITUTING* 

* PREVIOUSLY * 

* COMPUTED * 

* VALUES * 
*************** 



*************** 

* * 

* FIND MAXIMUM* 

* ELEMENT IN * 

* COLUMN * 

* * 
*************** 



ALL ROWS 
PROCESSED 



********* 



*** ** ********** 

EXIT 
*************** 



*************** 

* * 
* 

* SET KS TO 1 

* * 
*************** 



*************** 



*************** 



******************* 



*************** 



*************** 



50 V 55 

*************** 

* * *** 
♦DIVIDE ROW BY* * * 

* LEADING ****) BO* 

* COEFFICIENT * * * 

* * *** 
*************** 



SUftROUTlNC 6NAD0 



*************** 

FNTF.R 

*************** 



Hi************** 

* * 

* CALCULATE • 

* NUMBER OF * 

* ELEMENTS * 

ill * 

*************** 



*************** 



+ ADD MATRICES* 
+ R*A*B * 



*************** 



*************** 



SUBROUTINE GMSUB 



*************** 

ENTER 
*************** 



* CALCULATE * 

* NUMBER OF * 

* ELEMENTS * 
4 * 
*************** 



* SUBTRACT * 

* MATRICES * 

* R=A-B * 

* * 
*************** 



*************** 

EXIT 
*************** 



SUBROUTINE GHPRO 



*************** 



*************** 

FNTFR 
*************** 



** 



STEP POINTER* 

*g ***)« TO NEXT * 
* * * ELEMENT OF * 

** *OUTPUT MATRIX* 

*************** 



*************** 

* INITIALIZE * 

* POINTER AT 
♦FIRST ELEMENT 

* .IF OUTPUT 

* MATMX * 
*************** 



*YES 
* 
N ********* 



*************** 

* * 
* 

* SFT K'l 

* * 
*************** 



*************** 

* * *** 

* * * * 

* INCREMENT ****) A** 

* J BY I * * * 

* * *** 
*************** 



*******!*********** 



*************** 

** * * 

'* * * 

;i ^^***)* SET J=l 



********* 



** 



*#*****♦***+*** 



*************** 



A/ t ***)* CLFS^ OUTPUT 
* * LOCATION 



*************** 



*************** 

EXIT 

*************** 



******************* 



*************** 

* * 
* 

* SET 1=1 

* 

* * 
*************** 



*************** 



*** 



♦INCREASE K BY****1 A3* 

* 1 * * * 

* * *** 
*************** 



*************** 



*************** 



*************** 

* * 

* ACCUMULATE * 

* PROOUCT * 
*H(J,K) + BU,K)* 

* * A< J* I) * 
*************** 



*YES *** 

* * * 
I * M *| BO* 



*************** 



* INCHEMFNT I ****) A6* 

* W l * ***** 

* * *** 
*************** 



52 



SUBROUTINE GMTRA 



* * ** * +* * ** ** ** * 



**4*****#****** 



V 

*************** 

* INITIALIZE * 

* POINTER AT 

* FIRST ELEMENT 

* OF OUTPUT 

* MATRIX * 
*************** 



*************** 



*************** 



*************** 



* *NO **+ 

* * * 
LAST ELEHENT *) A2* 

* * * 

* * *** 
* * 

* * 
* YES 



EXIT 

*************** 



SUBROUTINE GTPftD 



*************** 



*************** 



*************** 

* INITIALIZE * 

* POINTER AT 
♦FIRST ELEMENT 

* OF OUTPUT 

* MATRIX * 
*************** 



*************** 



HA?**+|* GET ELEMENTS* 

* * * FROM A ANO B* 



*************** 



*************** 



* MULTIPLY * 

* TRANSPOSE OF* 

* A BY B * 

* * 
*************** 



* *NO 

* 
LAST ELEMENT ********* 



*************** 

EXIT 
*************** 



******************* 



*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ****! A2* 

* ELEMENT OF * * * 
•OUTPUT MATRIX* *** 
*************** 



SUBROUTINE MAOO 



*************** 
*************** ** J SET AMAY * 

* « * * POINTER AT 
ENTFR * *BO***l*FIRST ELEMENT 

* » » * FOR BOTH 

„♦«».,*,**,** » :,.j»j5jjii,„» 



* *NO 
BOTH MATR(CES * 

ST0RA6F M M0nE * """"**""*'* KEMENTS 

* + * ¥* # 

* * * ********#♦***** 



* AOD MATRIX A* 

* TO MATRIX B * 

* * 
**:***« ********* 



* CALCULATE * 

* RII.J) - * 

*a(1,j3+bii«j>* 



TO MATRIX B * J SWaN?' M S r" 

* *************** 



* • *NO 
PXIT * I * * LAST ELEMENT **•»***** 

* * * u. „ * * 

*************** * * * I 

* * . T 

* * * * 

* « YES » 

* * I 

* * * 

* * I 

******************* * * 

I * I 

* * 

V * 

*************** * 

CALCULATE * I cvit * * 

MTEST - * "'"" * * 



MS»*MSR * 
* * 

*************** 



*************** 



*************** 



******************* 
* 
* 
V 
*************** 

*** 



* SET OUTPUT * * STEP POINTER* 

* STORAGE MODE * J I9 s K§ X I nI1 ***** Bl J 
J MSR '° ,* IbOT^SXtRiSesJ *** 
*************** *************** 



MTFST grfater ********* 

THAN * * 



* SET MSR = 1 

* 

* * 
*************** 



* (*******#********* 



MTEST GREATER 
THAN 2 



* SET MSR =2 *» BO* 

* * * * 

* * *«* 55 
*************** 



SUBROUTINE MSUB 



*************** 

F^JTFR 
*************** 
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*************** 

** * SET ARRAY * 

* * * POINTER AT 
*RO*+*»*FIRST ELEMENT 

* * * FOR BOTH 

** * HATRKES ■" 

*************** 



* *NO 

Hi)TH MATRICES * 
OF SAME. 
STORAGE MODE * 



***** ********** 



**********Bl***j* 



GET BOTH 
ELEMENTS 



*************** 



* SUBTRACT * 
#MATRIX B FROM* 

* MATRIX A * 

*************** 



* CALCULATE * 

* R(I*J) * * 
*AI l.JJ-BUtJ)* 

* FOR ANY MSR * 
*************** 



*************** 
EXIT 

A* ******* ****** 



LAST ELEMENT 



********* 



******* ****** ****** 



*************** 

* * 

* CALCULATE * 

* MTCST - * 

* MSA+MSR * 
a * 
*************** 



*************** 



*************** 



*********:********** 



*************** 

* * 

* SET OUTPUT 

* STORAGE MODE 

* MSft=0 

* * 
*************** 



*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT *#**) Bl* 

* ELEMENT FOR * * * 
*BOTH MATRICES* *** 
*************** 



********* 



* SET MSK = 1 



*************** 



*( ***************** 
* 

* 35 

+ 

♦NO *** 
* * * 

*) BO* 
* * 



************ *** 

* * *** 

* * * * 

* SET MSR =2 *) BO* 

* * * * 

* * *** 
*************** 



SUBROUTINE NPRD 



4* ****** ******* 



*************** 
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*************** 

>: A «w TE : 

* R<JfK* + * 
*&U»K)4AU,(>* 

*************** 



********* 



********* 



*************** 

* * 

* SPECIAL CASE* 

* FOR DIAGONAL* 
♦BY DIAGONAL .* 

* R - A*B * 
** ** ******** *** 



******* #******#^< 



♦INCREASE I BY****I A9* 

* 1 * * * 

* * *** 

*************** 



******************* 



*************** 



*************** 



♦OUTPUT MATRIX* 

*************** 



******************* 



*************** 

* INITIALIZE * 

* POINTER AT 
♦FIRST ELEMENT 

* OF OUTPUT 

* MATRIX * 
*************** 



********* 



*************** 

* * 
* 

* SET K = I 

* 

* * 

*************** 



*************** 



♦INCREASE J BY****> AT* 

* 1 * * * 

* * *** 

*************** 



******************* 



*************** 

** * * 

* * * 
*A6***J* SET J = 1 

* * * 

** * * 

*************** 



********* 



*************** 
** * * 

* * * CLEAR OUTPUT 
#A7**+)* ELEMENT 

* * * LOCATION 

** * * 

*************** 



*************** 



*************** 



******************* 



******** **** *** 



******* ******** 



*************** 



♦INCREASE K BY****) A6* 

* i * * * 

* * *** 

*************** 



*************** 



*A9***)* GET ELEMENTS****) BO* 

* * * FROM A ANO B* * * 

** * * *** 

************ *** 



SUBROUTINE HT*A 



*************** 

FNTER 
*************** 



STORAGE MODE ********* 
GENERAL » * 



10 V 

*************** 

*MCPY * 

******* ******* 
* COPY MATRIX * 



*************** 



EXIT 
*************** 



******************* 



20 V 
*************** 

* INITIALIZE * 

* POINTER AT 
♦FIRST ELEMENT 

* OF OUTPUT 

* MATRIX * 
*************** 



*************** 



**;***)* GET ELEMENT * 

* * * FROM A * 

** * * 

*************** 



*************** 



* MOVE ELEMENT* 

* TO R * 



*************** 



* *NO «** 

* * * 
LAST ELEMENT *1 *5* 



EXIT 
*************** 



M 



*************** 



** ************* 



SUBROUTINE TPRO 



* *NO 

* 
LAST ELEHENT ********* 



DIAGONAL BY 
DIAGONAL 



********* 



*************** 
EXIT 

*************** 



******************* 



* CALCULATE » 

♦OUTPUT MATRIX* 

* FOR 01 AGONAL* 

* BV DIAGONAL * 
*************** 



*************** 

* * *** 

* STEP POINTER* * . * 

* TO NEXT ****> A6* 

* ELEHENT « * * 

* * *** 
*************** 



*************** 



*************** 



******************* 



*************** 

* * 

* INITIALIZE 
♦ARRAY POINTER 

* AT FIRST 

« ELEHENT * 
*************** 



*************** 

* * 

* CLEAR OUTPUT 

* ELEMENT 

* LOCATION 

* * 
*************** 



*************** 



*A6***I* GET ELEHENT * 



*************** 
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*************** 

* * *** 

* MULTIPLY * * * 

* TRANSPOSE OF****) BO* 

* A BY B * * * 

* * *** 
*************** 



SUBROUTINE HAT* 



*************** 

FNTFR 
*************** 



* * * LOCATE * 
*A1***I* FLFMENT IN * 

* * vtlUTPUT MATRIX* 
** * * 

*************** 



********* 



* CLEAR OUTPUT 

* LOCATION 



*************** 



*************** 

* * * 

* * GFT ELEMENT * 
<****)* JN A AND ITS* 

* * TRANSPOSE * 

* * * 
*************** 



*************** 



*ND *** 

* * * 

*» A4* 
* * 

*** 



*{***************** 

* 

* * 

* *NO *** 
LAST ELEMENT + * * 

IN OUTPUT *) Al* 
MATRIX * * * 

* * *** 



EXIT 
** ** **** * * ***** 



SUBROUTINE SAOD 



*************** 

ENTER 
*************** 



*************** 



* COMPUTE * 

♦VECTOR LENGTH* 
$ * 

*************** 



♦ADD SCALAR TO* 

* EACH ELEMENT* 

* * 
*************** 



EXIT 

*************** 
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SUBROUTINE SSUB 



*************** 
ENTFR 

*************** 



« COMPUTF * 
♦VECTOR LENGTH* 
* * 

*************** 



*************** 

* * 

* SUBTRACT * 
» SCALAR FROM * 

* EACH ELEMENT* 

* * 
*************** 



EXIT 

*************** 
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SUBROUTINE SNPV 



*************** 

ENTER 
*************** 



* COMPUTE * 
*VECTOR LENGTH* 



*************** 



"■MULTIPLY EACH* 

* ELEMENT BY * 

* SCALAR * 

* * 
*************** 



*************** 
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SUBROUTINE SDIV 



*************** 

FNTFR 
*•** *********** 



* COMPUTE * 
•VECTOR LENGTH* 



*************** 



********* 



*************** 



* SET VECTOR * 

* LENGTH TO 1 * 



*************** 



*(****•************ 



* DIVIDE * 

* ELEMENTS 3Y * 

* SCALA* * 

* * 
*************** 



EXIT 
*************** 



SUBROUTINE MOP 



*************** 



*************** 



*************** 

* * 

* INITIALIZE 

* POINTER AT 
•FIRST ELEMENT 

* OF ROW « 
*************** 



* * * GET ELEMENT » 

*A?»**I* FROM INPUT * 

* * * MATRIX * 
** * * 

*************** 



* ADD ELEMENT * 

* TO OUTPUT * 

* MATRIX * 

* * 
*************** 



LAST ELEMENT ********* 
OF ROW * * 



V 
*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ****> A2* 

* ELEMENT OF * * * 

* ROW * *** 
************ *** 



******************* 



EXIT 
*************** 



SUBROUTINE CADD 



*************** 

FNTFR 
4************** 



* INITIALIZE * 

* POINTER AT » 
♦FIRST ELFMENT * 

* OF COLUMN * 
*************** 



• * * GET ELEMENT * 
*AZ***1* FRtIN INPUT * 

• * * MATRIX * 

** * * 

*************** 



* ADD ELEMENT * 

* TO OUTPUT * 

* MATRIX * 

* * 
*************** 



LAST ELEMENT ********* 

OF COLUMN * * 



*************** 

* * **• 

* STEP POINTER* * * 

* TO NEXT ****) A?* 

* ELEMENT DF * * * 

* COLUMN * *** 
*************** 



******************* 



*************** 



SUBROUTINE SRMA 



4**4*********** 



*************** 



**** *********** 

* * 

* 

* SET J * 1 

* * 
*************** 



*************** 



*A2***>* GET ELEMENT5* 

* * * IN aOTH ROWS* 

** * * 

*************** 



********* 



*************** * 

* MULTIPLY * * 

* ELEMENT BY ************ 

* CONSTANT * ** 

* * ** 
*************** ** 



******************* 



I V 

*************** 

* MULTIPLY * 

* ELEMENT BY * 

* CONSTANT ANO* 

* ADO TU OTHER* 

* ROW * 

*************** 



*(**************** 



********* 



*************** 



******************* 



*INCREASE J 0Y****I AZ* 

* 1 * * * 

* * *** 

*************** 



SUBROUTINE SCHA 



*«#******««**** 

ENTER 
*************** 



*************** 



*************** 
** + * 

* * * GET ELEMENTS* 
sft?***l* IN I3HTH * 

* * * COLUMNS * 

** .* * 

**************# 



********* 



* MULTIPLY * * 

* ELEMENT BY ************ 

* CONSTANT * ** 

* * +* 

*************** ** 



******************* 



1 V 

4************** 

* MULTTPLY * 

* ELEMENT «Y * 

* CONSTANT AND* 

* 400 TO OTHER* 

* COLUMN * 
4************** 



* I**************** 



********* 



««*****«******* 



*************** 



******************* 



«*** ****** ** +*9 



♦INCREASE I BY****) A2* 

* 1 * * * 

* * *** 

*************** 



SUBROUTINE HINT 



*************** 

ENTER 
*************** 



* INITIALIZE 

* POINTER AT 
♦FIRST ELEMENT 

* qf sow * 

*************** 



CA?***)* GfT ELEMENT * 

.. * * IN BOTH ROWS* 

** * * 

*************** 



* INTERCHANGE * 

* FLEHENTS * 

* * 
*************** 



LAST ELEMENT ********* 
OF ROW * * 

* * * 



V 

*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ****! 42* 

* ELEMENT OF * * * 

* ROW * *** 
*************** 



******************* 



EXIT 

*************** 



SUBROUTINE CINT 



*****#*****++** 

FNTH* 
********* ****** 



**«**:********** 

* * 

* fNlTTUTZF 

* PQINTTK ftT 
*1 OUST CLFMFNT 

* UF C'JLUMN * 
*************** 



;= * * GET ELFMFNTS* 

.*•>**-*] * IN Bm H * 

* * COLUMNS * 



*************** 



*************** 



* IMTF'tCHAMGF * 

* tLE;HlrNTS * 

* * 
*************** 



l"(.f MLNT ********* 

;i"IUIMN * * 



*************** 

* * *** 

* STFP PHINTFP* * * 

* VI NSXT **+*> A2* 

* FlFMENT HF * * * 

* COLUMN * *** 
*************** 



******************* 



*************** 
FXIT 

*************** 



SUBROUTINE RSUM 



FNTPR 



*************** 



**«$**##«*#* «*# 



* SPT J=l 

* + 

********.*#*#*** 



*£ # * 

* * * GET ELEMENT * 

*A")***1* IJ IN INPUT * 

r * * MATRIX * 

** 4 * 



A * 

* ACCUMULATE * 
*SUM IN OUTPUT* 

* VICTOR * 

* * 

<:***#** £4 ****** 



A 



**+****+******* 





* *** 




* # * 


INCUFASE 


*#**) A3* 


J «Y 1 


* * * 




* *** 



****J|:#***##*ft4A 



******* £**£*#****#$ 
* 

V 
# 

*YES 
# 



* rNCREASn ****> A3* 

* I (IV 1 * * * 

* * *** 



******************* 



exit 

************ *** 



SUBROUTINE CSUH 



*************** 

ENTER 

*************** 



*************** 



*************** 



*************** 



*************** 



* * * GET ELEMENT * 
(rA3***|* IJ ]N INPUT * 

* * * MATRIX * 
** * * 

*************** 



* ACCUMULATE * 
*SUM IN OUTPUT* 

* VECTDR * 

* * 
*************** 



♦ YES 
* 
N ********* 



* INCREASE ****> A3* 

* I 8Y 1 * * * 

* * *** 
*************** 



******************* 



********* 



* INCREASE ****) A3* * 

* J BY 1 * * * * 

* * *** * 

*************** * 



******************* 



EXIT 
*************** 



SUBROUTINE RTA6 



*************** 

ENTER 
**** *********** 



*************** 



*************** 

* * 
* 

* SET I = 1 

* 

* * 

*************** 



B<i) outside ********* 

RANGE * * 



* ADD I TO ************ 

* SlL+l) * ** 

* * ** 
*************** ** 



******************* 



*************** 
*HA(1D * 

******* ******* 

* ADD ROW I OF* 

* 4 TO R * 



*************** 



ADO I TO S 
VECTOR 



*************** 



* < **************** 



********* 



*************** 

EXIT 

*************** 



******************* 



♦INCREASE I BY****j A3* 

* t * * * 

* * *** 
*************** 
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SUBROUTINE CTAB 



*************** 
ENTEft 

******** ******* 



*************** 



*************** 

* * 
* 

* SET I = I 
* 

* * 
*************** 



********* 



* ADD 1 TG ************ 

* sa+n * ** 

* * ** 
*************** ** 



******************* 



v 
*************** 

♦C40U * 

******* ******* 

* Aon cm. i of* 

* A TO ft * 

* * 
*************** 



*************** 



* [«******•******** 



********* 



*************** 



*************** 



******************* 



*************** 

* * *** 

* * * * 

♦INCREASE I BY****) A3* 

* l * * * 

* * *** 
*************** 



SUBROUTINE RSRT 



********** 4**** ** 

* * * 

FNTFK * *60*l 

* * + 

*************** ** 



* * 


20 


* * 




* *NQ 


*** 


* 


* « 


SORT 


*» A2* 


INDICATOR - * 


* * 



* BUILD « « * * * 

* ORIGINAL * *B1***I* GET SEQUENCE* 
*SEUUENCE LIST* * * * LIST NUMBER * 

* * ** * * 
*************** *************** 



* GET ROW IN * 

* INPUT MATRIX* 
*A CORRESP. TO* 

* * *SEQ. LIST NO.* 

*************** *************** 



« * * MOVE ROM TO * 

* SET 1-2 * *OUTPUT MATRIX* 

* * * R * 

* * * * 
*************** *************** 



* *YES * *N0 *** 

IK II GREATER * * * * * 

THAN OR EQUAL ********* * LAST ROW *) Bl* 

in 811-11 » * * * * * 



*************** 



* *************** 

*************** * 



* INTERCHANGE * 

* TWO ELEMENTS* 

* IN SORTING * 

* KEY VECTOR * 
4************** 



*************** 



* IN ORIGINAL * 
♦SEQUENCE LIST* 
*************** 



♦ I ***************** 

♦ 

V * 

*YES *** 

* * * 
= N *» BO* 

* * * 

* *** 



+INCREASF I BV****» A4* 

* 1 * * ♦ 

* * *** 

*************** 
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SUBROUTINE CSRT 



*************** 



ENTER * *BO*> SORT *1 A2* 

• * * * INDICATOR » * * * 

*************** », J.U1V.IU1, u ^^ 



*************** 



* BUILD * * * * • 

* ORIGINAL * *BL***I* GFT SEQUENCE* 
*SFOUENCE LIST* • * * LfST NUMBER* 

* * ** * * 
*************** *************** 



*************** 

* * * * 

* * * GET COL IN * 
'*. CLEAR SORT * « INPUT MATRIX* 

* INDICATOR * *A CORRESP. TO* 

* * *SEO- LIST NO.* 
*************** ************»}* 



*************** 

* * *ccpy * 

* . _ * ******* ******* 

* SET J - 2 * « MOVE COL TO * 

* * * OUTPUT MAT R* 

* * * * 
*************** *************** 



S * «NO *** 

********* * LAST COLUMN *» BJ* 



*************** 



*************** 

* * 

* INTERCHANGE * 

* TWO ELEMFNTS* 

* IN SORTING * 

* KEY VECTOR * 
*************** 



*************** 

* * 

* INTERCHANGE * 

* TWO ELEMENTS* 

* IN ORIGINAL * 
♦SEQUENCE LIST* 
*************** 



*************** 



* ( ***************** 



♦YES *** 



♦INCREASE I BY****I A4* 
* 1 * * * 

+ * *** 

*************** 



SUBROUTINE RCUT 



*************** 

ENTER 
*************** 



*************** 

* * 

* INITIALIZE 
*ARRAY POINTER 

* AT FIRST 

* ELRIFNT * 
*************** 



#42***)* get ELEMENT * 
* * * * 

** * * 

*************** 



* *YES 

ELEMENT IN * 
COLUMN L OR ********* 

BFLOW L * * 



*************** 



HOVE ELEMENT************ 
TO MATRIX R * ** 



*************** 



******************* 



*************** 



*************** 



*(**************** 



LAST ELEMENT ********* 



*************** 



******************* 



*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT **#*( A2* 

* ELEMENT * * * 

* * *** 
*************** 



SUBROUTINE CCUT 



*************** 

ENTER 

*************** 



*************** 



AT FIRST 

* ELEMENT * 

*************** 



*A2***»* GET ELEMENT * 

* * * * 

** # # 

*************** 



40 * * 
* * 

♦ YES 



PLEMFNT IN 
COLUMN L OR ********* 

SIGHT OF L * + 



*************** 



*************** 



«**********#******* 



*************** 



* MUVE ELEMENT* 

* TO MATRIX S * 



*************** 



*(**************** 



LAST ELEMENT ********* 



*************** 



******************* 



V 
*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ***** A2* 

* ELEMENT * * * 

* * *** 
*************** 



SUBROUTINE RUE 



******** ******* ** ********* ****** 

* * * * 

FNTFR * *BO*)* EXIT 

* * * * 
*************** ** *************** 



* * ** * # 

* SfT STORAGF * * * *RESET STORAGE 

* MODE ANU * *R1***)* MODE AND. 

* NUMBER OF * * * * NUMBER OF 

* ROWS FHR A * ** * ROWS FOR A * 
*************** *************** 



V V 

*************** *************** 

** * INITIALIZE * * * *** 

* * *Al<KAY POINTER * * * * * 
*4 ? ***>* AT FIRST * * GET ANOTHER *+** ) A2* 

* * * ELEMENT OF * * COLUMN * * * 
*« * COLUMN * * * *** 

*************** *************** 



*************** 

** * * 

n * * * 

* A }***)* GET ELEMENT * 

* * * * 

ft* * * 

*************** 



*************** 



*************** 



LAST ELEMENT ********* 

"IF COLUMN * * 



*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ****» A3* 

* ELEMENT * * * 

* * *** 
*************** 



******************* 



*YES 
BOTH INPUT * 

MATRICES ********* 

PROCESSED * * 



^RFSET STORAGE 

* MODE AND 

* NUMBER OF 

* RUHS FOR B * 
*************** 



******************* 



* *NO *** 

LAST COLUMN *) Bl* 



* *** 

* YES * * 

********** ) * BO* 



SUBROUTINE CTIf 



*************** 



*************** 



$************** 

* * 

* SFT STORAGF 

* MODE AND 

* NUMBER OF 
♦COLUMNS FF1R A * 
*************** 



+* 



*************** 



INITIALIZE 
*A2***I*AR«AY PUINTFR 
* * * AT FIRST 
** * ELEMENT * 

*************** 



*************** 



*/!n***)* GET ELEMENT * 



*************** 



*************** 



* MOVE ELEMENT* 

* TO MATRIX R * 



*************** 



LAST ELEMENT 



********* 



*************** 

* * ** + 

* STEP POINTER* * * 

* TO NEXT ****) A3* 

* ELEMFNT * * * 
4 * *** 
*************** 



******************* 



BOTH INPUT 

MATRICES 
PROCESSED 



********* 



*************** 



*************** 



******************* 



*************** 

* * 
*RESET STORAGE 

* MODE AND 

* NUMBER OF 
♦COLUMNS FOR B * 
*************** 



SUBROUTINE KCPY 



*************** 



*************** 



*************** 



* COMPUTE * 
♦VECTOR LENGTH* 



*************** 



* COPY MATRIX * 
« ELEMENT BY * 

* ELEMENT * 

* * 
** ************* 



*************** 



SUBROUTINE XCPY 



*************** 
ENT6K 

*************** 



* COMPUTE * 

* LIMITS L2.K2* 



*************** 



*************** 

* * 

* INITIALIZE 

*ARRAY POINTER 

* AT L,K IN 

* MATRIX A * 
*************** 



*A-1***)* GET ELEMENT * 

* * * * 

** * * 

*************** 



* HOVE ELEMENT* 

* TO MATRIX R * 



*************** 



LAST ELEMENT ********* 
(L2.K2) * * 



*************** 



******************* 



*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ****» A3* 

* ELEMENT * * * 

* * *** 
*************** 



SUBROUTINE RCPY 



************+*# 
ENTER 

*************** 



* INITIALUE 

* POINTER AT 
*F1RST ELEMENT 

* OF ROW * 

*************** 



» * * GET ELEMENT * 
*A?***)* FROM INPUT * 
* * * MATRIX * 



*************** 



* MOVE ELEMENT* 

* TO OUTPUT * 

* VECTOR * 

* * 
*************** 



LAST ELEMENT ********* 
OF ROW * * 



*************** * 

* * *** * 

* STEP POINTER* * * * 

* TO NEXT ****] A2* * 

* ELEMENT OF * * * * 

* ROW * *** * 
*************** * 



******************* 



*************** 



*************** 
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SUBROUTINE CC»Y 



*************** 

FNTFR 
*************** 



*************** 



*FIRST ELEMENT 

* nF COLUMN * 
*************** 



ft * * GET ELEMFNT * 
*4?***)« frqh INPUT * 
* * * MATRIX * 

** # * 

*************** 



***#****#****** 

* * 

* MOVE ELEMENT* 

* TO OUTPUT * 

* VECTOR * 

* * 

**4* *********** 



*************** 

* * *** 

* STEP POINTER* * * 

* TO NFXT ****) A2* 

* ELEMENT OF * * * 

* COLUMN * *** 
*****«*#******* 



******************* 



EXIT 

*************** 



SUBROUTINE DCPY 



*************** 



*************** 



*************** 

* * 

* INITIALIZE 

* POINTFR AT 
*FIRST ELEMENT 

* * 
*************** 



> * * GET DIAGONAL* 
-Ay***l* ELEMENT FROM* 
* * * INPUT MATRIX* 

** * * 

** ********** *** 



* MOVE ELEMENT* 

* TO OUTPUT * 

* VECTOR * 

* * 
*************** 



LAST ELEMENT ********* 
OF DIAGONAL * * 



*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ****! A2* 

* ELEMENT OF * * * 

* DIAGONAL * *** 
*************** 



******************* 



*************** 

EXIT 
*************** 



SUBROUTINE SCLA 



*************** 



*************** 



*************** 



*************** 



* REPLACE BY * 

* SCALAR * 

* ELEMENT 8Y * 

* ELEMENT * 
*************** 



EXIT 

*************** 



SUBROUTINE OCLA 



FNTC-ft 



* INI riALl/.E 

* PQINTFtl AT 
*MHST fcLEMfcNT 

* * 
###*:###*** ***** 



fl?***l* [jFI III AGONAL* 

, * # iTtFMCMT * 

** * * 

*************** 



* * 

* KEPLACF * 
v 01 AGONAL * 

* ' LFMFNT ril TH* 

* SCALES * 
*************** 



********>> 



**#****#*#*#*** 

* * *** 

* STFP POINTER* * * 

* TO NPXT ****) A2* 

* Ft FMENT UF * * * 

* DIAGONAL * *** 
**** *********** 



a****************** 

* 
ft 

* 
V 

*******####**** 

* 
FXIT * 

* 
*************** 



SUBROUTINE HSFit 



*************** 

ENTER 
*************** 



*************** 

* * 

* INITIALIZE 

* ARRAY POINTFR 

* AT FIRST 

* ELEMENT « 
*************** 



*A2***I* GET ELEMENT * 



*************** 



* *YES 

* 
OUTPUT MATRIX ********* 
GENERAL * * 



*YES * 

ELEMENT IN * * 

LOWER ********* 

TRIANGLE * ** 



*(***************** 



* *YES 

UFF-DIAG * 

ELEMENT OF * 
DIAGONAL * 
OUTPUT MATRIX 



*************** 

* * 

* FORM ELEMENT* 

* OF OUTPUT * 

* MATRIX * 

* * 
*************** 



* | ***************** 



LAST ELEMENT ********* 



*************** 



*************** 



******************* 



V 
*************** 

* * *** 

* STEP POINTER* * * 

* TO NEXT ****! A2* 

* ELEMENT * * * 

* * *** 
*************** 



SUBROUTINE MFUN 



*************** 



*************** 



*************** 



* COMPUTE * 
•VECTOR LENGTH* 



*************** 



*************** 

* APPLY * 

* FUNCTION TO * 

* MATRIX * 

* CLEMENT BY * 

* ELEMENT * 
*************** 



EXIT 
*************** 
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FUNCTION RECP 



*************** 

ENTER 
*************** 



* *VES 

* 
ZERO ELEMENT ********* 



* CALCULATE ************ 

* RECIPROCAL * ** 

* * ** 

*************** ** 



***** ************** 



* SET EQUAL TO* 

* INFINITY * 



*************** 



* ( **************** 



EXIT 

*************** 



SUBROUTINE LQC 



*************** 
ENTER 

*************** 



********* 



* COMPUTE * * 
♦SUBSCRIPT FOR************ 

* GENERAL * *+ 

* MATRIX * ** 
*************** ** 

** 
** 



******************* 



MS EQUAL TO 1 ********* 



* ** 

* ** 
V ** 

*************** ** 

* * *** ** 

* COMPUTF * * * ** 
♦SUBSCRIPT FOR****) A6* ** 

* SYMMETRIC * * * ** 

* MATRIX * *** ** 
*************** ** 

** 



******************* 



* COMPUTE * 
♦SUBSCRIPT FOR* 

* DIAGONAL * 

* MATRIX * 
*************** 



* ( **************** 



** *************** 

*A6*>* EXIT 

* * * 

** *************** 



SUBROUTINE ARRAY 



*************** 
FNTFR 



IS MODE - I ********* 



* CONVERT FROM* « 

* SINGLE TO ***»***«*»** 

* DOUBLE * ** 

* DIMENSION * ** 
*****+♦+****#** ** 



******************* 



* CUNVERT FROM* 

* DOUBLE TO * 

* SINGLE * 

* DIMENSION * 

*************** 



+ 1 **************** 



*************** 



*************** 



SUBROUTINE QUADR 



*************** 

******** ******* ** * « 

* * * * * 
ENTER * *BO***»*ADD RESULT TO* 

* * * * SUM * 
*************** ** * * 

* *************** 



*************** 

* * * * 

* * * STEP TO NEXT* 
*SFT SUM EQUAL* * POINT IN + 

* Tno* * FUNCTION * 

* ♦ * TABLE * 
*************** *************** 



* *NO ** 

ARE THERE * * * 

LESS THAN TWO *****+****R2* ) ARE THERE * ********* 

FUNCTIONAL * ** * " — 

* VALUES * * ** 



* *************** 

* * * * * 

* * * * USE E0.(2J * 

* SET ERROfl * * * AND ADD TO * 

* CODE TO 1 * * * SUM * 

* * * * * 
*************** * *************** 



*************** 



*************** * *************** 



******************* ******************* 

* * 

V V 

^1 * * * * 

* * * * 

* DOES THE *NO * *NO 
SPACING OF * * * 

FUNCTION ********* * ARE THERE 3 ********* 
VALUES EOUAL * * * POINTS LEFT * * 

* ZERO * * * * * 



*************** * 

* * * * * * 

* * * * USE EQ.(3) * * 

* SET ERROR * * * and ADD TO ************ 

* CODE TO 2 * * * SUM * ** 

* * # * * ** 
*************** * *************** ** 



* * ******************* 

* * * 

* * V * 

* * * * * 
v * * * * 

*************** * * *no * 

* * * * y 
EXIT * * * ARE THERE 2 ********* 

* * * POINTS LEFT * * 
*************** * * * * 



******************* 



* *YES *** * * 
ARE THERE * * * * USE E0.<4) * 
LESS THAN *) B2* * AND ADD TO * 

FIVE POINTS * * * * SUM * 

* REMAINING * *** * * 
* * *************** 



*{ ***************** 



39 V * 

*************** y 

* * *** *************** 

* USE 5 POINT * * * * 

* FORMULA ****) BO* * EXIT 

* EQ. (1 > * * * * 

* * *** *************** 
*************** 



SUBROUTINE SHPSN 



*************** 



*************** 



*************** 



♦ * 



nivine * 

♦BO***»*INTEKVAL INTO* 

* * * TWO * 
»* * SUBINTERVALS* 

*************** 



*************** 



*************** 



# *ND 

* AUE THE UPPER * 

AND LOhFR 

* INTEGRATION * 
♦LIMITS F'JUAL 



********* 



*************** 



PRESET 
CONSTANTS 



*************** 



*************** 



SFT ERROR 
CODE TJ I 



*************** 



*************** 

* CALC VALUE * 

* FUR INTEGRAL* 

* USING EQ.I1>* 

* FOR Z * 

* INTERVALS * 
*************** 



*************** 



*************** 



* *************** 

* ** # * 
** * * DIVIDE * 
**B3***)* INTERVAL BY * 
** * * IWO * 

* ** * * 

* *************** 



************ ******* 



IS THE GIVEN 

ACCURACY 

LIMIT ZFRd OR 

* NEGATIVE * 



********* 



*************** 

* * 

* COMPUTE NEW * 

* VALUE FOR * 

* INTEGRAL * 

* USING EQ.UI* 
*************** 



*************** 



*************** 



n * * 
* * 

♦NEW VALUE -+YES 
OLD VALUE 
LESS THAN 
REUUlKtD 

♦ TOLERANCF 



********* 



*************** 



*************** 



HAVE THE *NO 
MAXIMUM * 

NUMBER CJF 
ITERATIONS ■ 
BFFN OUNE * 



******************* 



* IS THE GIVFNNC1 

MAXIMUM ♦ 
ITERATION 
COUNT LFSS * 

* THAN I * 



*************** 



SET ERROR 
CODE TQ ** 



************ 



*************** 



******************* 



*************** 



*************** 



* SET ERROR * 

* CODE TO * 

* * 
*************** 



*( **************** 



*************** 

EXIT 
*************** 



♦ COMPUTE * 

♦ NUMBER OF ** 

* INTERVALS ♦ 

* ♦ 
*************** 



SUBROUTINE RK1 



*************** 

ENTER 
*************** 



•NO 



»BO*» 



XN1 * 
XI FINAL) 



* YES 

* 
* 
* 



**********CO*l 



•ves 
is this yn ********* 



XIFINAL) LESS 

THAN m FOUAL 
TO XIINITIAL) 



********* 



30 V 

*************** 

* * *** 
•SET XIANSHERl* * * 

* . X(FINAL), ****> C9* 

* YIANSWER) - * * . * 

* YN1 * •** 
*************** 



100 V I 

*************** * 

* x,AN iS!^ * ..♦*....**.* 

* Y(ANSWER) = * ** 

* YNl * ** 

*************** ** 

** 
** 



* XI ANSWER) = * 

* Xt INITIAL) , * 

* YIANSWERI = * 

* V( INITIAL) * 

*************** 



******************* 



*Y6S 



XN1 LESS THAN 

Xt FINAL) 
* * 



******************* 



80 V * 

*************** * 



*************** 



*************** 



*************** 



40 V 

*************** 

* * 

* » 

* SET HINEWI ** 
*X(F1NALI-X<N1* 

* * 
*************** 



* *************** 

* ** * * 
** * * INTERPOLATE * 

**C3«**I* TO FIND NEW * 

** * * TIHE STEP * 

* ** * * 

* *************** 



******************* 



IS STHP SI7E 



********* 



*************** 

* * *** 

* INTEGRATE * * * 

* Y(ANS) - ****> C9* 

* YIXN+HI * *** 
*************** 



*************** 

* * 

* * 
♦INTEGRATE YN1» 

* =. Y(XN*H) » 

* * 
*************** 



SET ERROR 
CODE TO 1 



*************** 



******************* 



* *YES 

* 00 YNl AND YN ********* 

* SPAN YIFINALI * * 



YNl = 
Y( FINAL ) 



********* 



*************** 

* * 

* X(ANSWER) * * 

* X( INITIALI. * 
♦Y(ANSWER) - 0* 

* * 
*************** 
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**********•••*• 



* YN = YNl, XN* 

* - XN1 * 



*************** 



150** 

** 

*YES *** ** 

* * * ** 
♦ I CB* ** 

* * * ** 
* *** ** 

** 
** 



*************** 



*************** 



*************** 

* * 

* * 
♦INTEGRATE YNl* 

* . YIXN) * 

* * 
*************** 



*************** 



*** ** 



1+1 



****) C3* ** 
***** 
* * *** ** 

*************** ** 

** 



*** **************** 



*(***************** 



******************* 
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* XN = 

* X1INITIAL). 



* 150 V 

* ^ *************** 

' nnc<; yn =. * Y * S ***** **** * XIANSHERI - * 

^W " «*i C °* *"*"'* YCANs'SeR. - * 

. yifVal. ** *~ « ;**.,,..*...... 



* NO 



*( ***************** 



*************** 



*** 



♦INTEGRATE YNl****) 80* 
* = YIXN) * *„/ 

*************** 



I » 1 



*•**) C3* 



* * * 
*C9*>* 

* * * 



*************** 



*************** 
EXIT/ 

*************** 



SUBROUTINE RKZ 



*************** 

*«*************: ** * * 

■* * * * * 

FNTfR * *I10***I*VECT0RII) * Y* 

* * * * * 

*************** ** * * 

* *************** 



*************** * * 

* * * *no 

* X = * * * 

#X( INITIAL), Y* * I « N ********* 

« = Y(INITIAL)* * * * 

* * + * + 
*************** * * * 



*************** 

« * 

* * 

* t = I * * FXIT 

* * * 

* * *************** 
*************** 



******************* 



*************** *************** 

+ * * * *** 

* * * * * * 
1* J = 1 * * I ■ 1+1 ****) A3* 

* * * * * * 

* * * * *** 

4************** *************** 



9************** 



* ftlNTEGRATF Y ** 
A',***)* Y(X+H> U^INii* 

* * F.O . ( ? ) « t 3 I * 
**■ * * 

*************** 



* X = X + H * 

* * 

* * 
*************** 



*yF<i *** 

* * * 
*l BO* 

* * * 
* *** 



*************** 

* * *** 

* * * * 

* J = J*-l ****) A** 

* * * * 

* * *** 
*************** 



SUBROUTINE RK3 



*♦*+**** #*#***# 



********* 



********* ****** 

* * 

* SET INITIAL * 

* VALUF FflH * 

* XtY(I) * 

* i --> i , i M s * 

* * ** * ** * * *** ** * 



*************** 

* * *** 

* * * * 

* ADO 1 TO LL ****) A3* 

* * * * 

* # *** 
*************** 



******** ******* 

** * * 

s * * * 

>A3***1* JJ=1 * 

S= * # * 

*# * * 

*************** 



** 



***#******#**** 



CHMPJTP * 

v 44.***) *'UJMGF-KUTTA K* 

« * VALUES (F-O. * 

** * I-'t) * 

#*****«******** 



******************* 



*************** 

* * 

* * 

* LL~l * 

* * 

*************** 



*************** 

EXIT 

*************** 



*************** 



*************** 



onisir- k STSP5 ********* 



*************** * 

* ***** 

* * * * * 

* ADD I TO J J ****) AA* * 

* * * * * 
4 * *** * 
*************** * 



******************* 



*************** 

* * *** 

* STORE NEXT * * * 
*Sr-T HF VALUES****) BO* 

* IN VAl. * * * 

* * *** 

*************** 



SUBROUTINE FOR IP 



*************** 

ENfh-R 
**** **** ******* 



*************** 



*RO***)* SET C * 1, S* 
** * - 0, J = 1 * 

** * * 

*************** 



** *************** 
+ * * 

*C0*>* EXIT 
* * * 

** *************** 



*************** 



*************** 



*************** 

* * 

* SET U2 - 0, * 
*>* Ul = 0, I - * 

* 2*N * 

* * 
*************** 



+ * * COMPUTE NEW * 

*C1***I* COS. Q = * 

* * * Cl*C-Sl*S * 

** * * 

*************** 



* IS NO OF *N0 

POINTS * 
DEFINING 
INTERVAL IN) * 

* NFRATIVE * 



********* *R2*** 



*************** 

*EUN * 

******* ******* 

»* OBTAIN NEXT * 
+ VALUE OF FUN* 
* * 

*************** 



*************** 

* + 

* CUMPUTE NEW * 

* SIN. S = * 

* Cl*S+Sl*C * 

* * 
*************** 



*************** 



*************** 



*************** 



SET EAR03 
CODE * ?. 



************ **^t 



* COMPUTE UO * 

* EQUATION (31* 



*************** 



* REPLACE C BY* 

* t) * 

* * 

*************** 



*************** 

FXIT 
*************** 



*************** 

* * 
♦REPLACE U2 9Y* 

* Ul AND Ul df* 

* UO + 

* * 
*************** 



*************** 

* * 

* * 

* ser j = j+i *** 

* * 

* * 
*************** 



*** 
* * 
) .u* 



******************* 



********* 



*************** 



*************** 



*************** 



* RET EARQR * 

* CODE * I * 

* * 
*************** 



*** 

* * 
*1 B2* 

* * 
*** 



*************** 



*************** 
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*************** 

* COMPUTE * 

* FOURIER * 

* COEFFS - * 

* A(J), B(J). + 

* EQ. 14), (5) * 
*************** 



******************* 



*************** 



*************** 



*NC *** 

* * * 
M+l *) CI* 

* * * 
* *** 



*************** 



* COMPUTE CI ****! 80* 

* FOUATIQN ( 11* * * 

* + **+ 
*************** 



* MODIFY **+*) CO* 
•CONSTANT TERM* * * 

* * *** 
*************** 



SUBROUTINE FORIT 



*************** 

ENTFJft 
******** ******* 



* ************** 



*BO**+)* SET C » 

* * * « 0, J 



1* s* 

=- 1 * 



*************** 



*************** 



*********** **** 



*************** 



* SET EKftDR * 

* CODE * + 

* * 
*************** 



* * * SET 02 - 0, * 
*B1***J* Ul - 0, I - * 

* * * 2*N * 
** * * 

*************** 



* * * COMPUTE NEW * 

*ci***»* COS. Q = * 

* * * Ci*C-Si*S * 

** * * 

*************** 



* IS NO OF *NO 
POINTS * 
DEFINING 
INTERVAL <N) * 
f NFGATIVE * 



*************** 

** * + 

* * * OBTAIN I^TH * 

**********B2***|* VALUE OF * 

** * * FUNCTION * 

* ** * * 

* *************** 



*************** 

* * 

* COMPUTE NEW * 

* SIN. S - * 

* Cl*S+Sl*C * 

* * 
*************** 



*************** 



*************** 



*************** 



SET ERROR 
CODE = 2 



*************** 



* COMPUTE UO * 

* EQUATION (3)* 



*************** 



* REPLACE C BY* 



*************** 



*************** 



*************** 



*************** 



*************** 



*************** 

* * *** 

* SET J = J+l ****) Hi* 

* * * * 

* * *** 

*************** 



******************* 



*************** 



********* 



* SET I = I-l * 



** *** ***** ***** 



*************** 



*************** 



*************** 



*************** 



60 V 

*************** 

* COMPUTE + 

* FOURIER * 

* COEFFS - * 

* MJ>, B(J>. * 

* EQ.<4).<5> * 
*************** 



******************* 



COMPUTE SI * 
EQUATION (2)* 



*************** 



*************** 



* COMPUTE CI ***** BO* 

* EQUATION Ml* * * 

* * *** 

*************** 



* MODIFY ****) CO* 
♦CONSTANT TERM* * * 

* * *** 

*************** 



SUBROUTINE C*MN4 



***** ********** 

F-NTFH 
**«*#*****#**** 



50 * * 110 
* * 

** * *ves *** 

* * * * * * 

*B0*) X GRfcATER *l AT* 

* * * THAN 1 * * * 
** * * *** 



*************** 



* *yes *** 

* * * 

X EQUAL TO I *1 A9* 



*************** 



*************** 

* * 

* * 
+ SFT GX T'l I * 

* * 

* * 
*************** 



60 * * 

* * 

* *N0 
X WITHIN * 

0.000001 OF ********* 

BEING * * 

* NEGATIVE * * 



*NfJ * + * 

* * * 

*> BO* 



*************** 



************** f< 



*************** 
EXIT 

************* **: 



************** ;< 



*************** 



* MULTIPLY OX * 

* <1Y X * 

* * 

*************** 



******************* 

* 

V 

* 110 

* 

*YES *** 
* * * 
*> A7 + 
* * 
*** 



X GftFATES 
THAN I 



#yr-s *** 



BO V 

**************:& 

* « 

* f 

* OIVIliE GX 1Y* 

* X » 

* * 

*************** 



* * cniPUTF. GY * 

*A7***)+ USING * 

* * * EQUATION 4 * 

** * * 

*************** 



* * *** 

* * * * 

* 400 1 TO X ****) 85* 

* * * * 

* * *** 

*************** 



*************** 

* * 

* MULTIPLY GX * 

* f*Y GY TO ♦ 
♦11RTA1N ANSWER* 

* * 
*************** 



120 * 

V 
** *************** 

y * * 

*A<?*>* EXIT 

V * * 
** *************** 



SUBROUTINE LEGEN 



********* ****** *: 

* * * * * 

ENTER * *BO***l* SET PNX TO 0******#***** 

* * * * • * 
*************** ** * * * 

* *************** * 



* * * PRESET * 
*B1***)*COMPUTATIONAL* 

* * * VARIABLES * 
** * * 

*************** 



********* * SET J * l * 

* * * 

* * * 

* *************** 



* ******* PERFORM * 

* SET ERROR ****> flO* **B3***)* RECURRENCE * 

* cone to 1 * * * ** * * eo.il) for * 

* * *** * ** * pfj) * 

*************** * *************** 



********) * j * n-I ********* 



*************** * 

* * *** * * * *** ** 

* **** * ***** 

* S£T ANSWER ***** B7* * * ADO 1 TO J ****) B3* ** 

* P(N> TO X * * * * * ***** 

* * *** * * * *** ** 
*************** * *************** ** 



******************* ******************* 

* * * 

V * * 

17 * * 20 100 V * 

* * *************** * 

* *Y(5S *** * * * 
IS X BETWEEN * * * * * * 
MINUS 1 AND *» Bl* * SET ANSWER * * 

PLUS l * * * * P(NI TO P(J)* * 

* # *** * * * 

* * ***** ********** * 



* SET ERROR *+**) BO* *B7***J* SET ERROR * 

* COOE TO Z * * * * * * COOE TO * 

* * *** ** * * 
*************** *************** 



*{ ***************** 

* 

* 

150 * 

V 

*************** 

* 
EXIT * 
* 
*************** 



SUBROUTINE 6ESJ 



*************** ** * ^ ** * *YFS 

* * * * * * * * HAS UPPER * 

ENTER * *BO***l* SET ERROR * *CO*l LIMIT FOR H ********* 

* * * * CODE TO 4 * * * * BEEN REACHED * + 

*************** ** * * ** * * * 

* *•**••********, X * * * 



*************** 

* * *** 

* * * * 

* SET BJ TO * * EXIT * * STEP M ****) 85* 

* * * '* * * * * 

* * ************** * * * *** 
*************** *************** 



*********** *****.*** 



* * * * * * USE CURRENT * 
N NEGATIVE **********B2***1* SET ERROR * * VALUE OF BJ * 

* ** * * CODE TO * * FOR ANSWER * 

* * ** * 4 * * 

* * * *************** *************** 



* *************** *************** 

* * * * a * * 

* * * * COMPUTE 4 * * 

* SET ERROR * * * STARTING AND* * SET ERROR * 

* CODE TO 1 * * * UPPER LIMITS* * CODE TO 3 * 

* * * * fqr m f * * 
*************** * *************** *************** 



200 * 
V 
** *************** 

* * * * * * * 

EXIT * * * INITIALIZE M* *C**I* EXIT 

* * * * * * * 
*************** * * * ** *************** 

* *************** 



******************* 



X NEGATIVE OR **********B5***1* SET STARTING* 

ZERO * ** * * VALUE FOR F * 

* * * ** * § 

* * * *************** 



*************** * 

* * * * * 

* * * * COMPUTE * 

* SET ERROR * * * VALUES FOR * 

* CODE TO ?. * * * ALL F * 

* * * * * 
*************** * *************** 



* *************** 

* * * 

* * *COMPUTE ALPHA* 
EXIT * * * BY SUMMING * 

* * *F«S (EQUATION"" 
*************** * * 3) * 

* *************** 



******************* 



* COMPUTE NEXT* 

* VALUE FOR J * 

* 8ESSEL * 

* FUNCTION * 

*************** *************** 



40 * * 200 

* * 

*VES *** * NEXT VALUE*VES *** 

* * * * WITHIN * * * 

* TOLERANCE OF *> C4* 
* PREVIOUS * * * 
* VALUE * *** 

• * 
* * *** 

* NO * * 
*******#**)* CO* 



N 


LESS THAN 




*1 B2* 




NTEST 


* 


* * 




4 




*** 


* * 








* * 




*** 




* NO 




* * 




********** )* BO* 








* * 








*** 
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SUBROUTINE BE5V 



*************** 



******* ******** 



*************** 



*B0***)# ADO TERM TO * 
♦ * * SUM OF TERMS* 



***** ********** 



********* 



* *NO 

* 
N LESS THAN 



********* 



* *YES 

DONE 16 TERMS 



********* 



********* 



*************** 



*************** 

* * *** 

•multiply term* * * 

* to produce ****! bo* 

* next term of* * * 

* eq.ui * *** 

*************** 



* «* 

* ** 

110 V ** 

*************** ** 

* * *** ** 

* SET ANSWER * * * ** 

* EQUAL TO ****) C9* ** 

* SECOND 1"EKM * * * ** 

* * *** ** 
*************** ** 

** 



******************* 



******************* 



*************** 



***# **** ******* 



*************** 

* * 
♦COMPUTE FIRST* 

* TERM OF * 
♦EQUATION UOI* 

* * 
*************** 



* SET ANSWER * * 

* EQUAL TU ************ 

* FIRST TEKM * ** 

* * ** 
«** ***** ******* ** 



******************* 



******************* 



SET ERROR 

CORE TO 



*************** 



*************** 



*F)4***|* ADD TERM TO * 
* * * SUM OF TERMS* 



*************** 



*K = l v Y(A) =* 

* Y(0>. Y(BI = ♦ 

* Yd) * 

* + 
*************** 



X LESS THAN 



********* 



* *YES *** 

DONE 16 TERMS *) CO* 



* COMPUTE * 
)* RECUKRENCE * 

* RELATION * 

* E3.ll) * 
*************** 



*************** 



*************** 

* * *** 
♦MULTIPLY TERM* * * 

* TO PRODUCE ****) B4* 

* NEXT TERM OF* * * 

* EQ* UOI * *** 

*************** 



********* 



*************** 



*************** 



* 30 

* *************** 

* ** * COMPUTE * 
** * * COEFFICIENTS* 
**B7***I* FUR SINt COS* 
** * * TERMS * 

* ** * EQ,(4)-I7) * 

* *************** 



* ** 

* ** 
V 140** 

*************** ** 

* * *** ** 

* STEP INDEX * * * ** 
♦FOR NEXT TERM****) C5* ** 

* OF EQ. (II * * * ** 

* * *** ** 
*************** ** 

** 
** 



******************* 



******************* 



*YES *** 

* * * 
*) B7* 

* * * 
* *** 



*************** 
♦COMPUTE FIRST* 

* TERM OF * 

* RECURRENCE * 

* RELATION * 

* EC. 12) * 
*************** 



•SET ANSWER TO* 

* CURKENT TERM* 
♦OF RECURRENCE* 

* RELATION * 
*************** 



*( **************** 



40 V 

*************** 

* * *** 

* COMPUTE * * * 

* INITIAL TERM****) 60* 

* OF EQUATION * * * 

* m * *** 

*************** 



*************** 

* COMPUTE * *** 

* SECOND TERM * * * 
*0F RECURRENCE****) CO* 

* RELATION * * * 

* EU,<3> * *** 
*************** 



*************** 



*************** 



SUBROUTINE BE5 I 



*************** 



*************** 



*************** 
*+ * * 

* * * COMPUTE NEXT* 
•80***1* TERM OF * 

* * * EO.(L) * 
** * * 

************** * 



*************** 



*************** 



*************** 



********* 



*************** 



*Cl***>* COEFFICIENT * 

* * * TO COMPLETE * 

** * EQ.I2I * 

*************** 



* SFT FRROR **** » C2* 

* CODE TO 1 * * * 
w * *** 
************** * 



* * * 
**********C2* 1 * 

* * * 



*************** 



*************** 



******************* 



*************** 



********* 



****! A9* 



*************** 



* SET ERROR ***+) C2* 

* CUOE TO 2 * * * 

* * *** 
*************** 



** * +COMPUTE FIRST* 
**34***>* term OF * 
** * * £0.(21 * 

• ** * * 

* *************** 



******************* 



*************** 

* * 
♦SET TOLERANCE* 

* = .000001, * 

* SFT ERROR * 

* CODE TO * 
*************** 



*************** 



- TERM, 
« 1 



*************** 



* *YFS **» ** 

* # * * * 

IS X GREATER *) B<t* *B6*> 

THAN 12 AND N * * * * * 

* * *** ** 



* *VES *** 
IS DIFFERENCE * * * 
BETWEEN TERMS *) CI* 

LESS THAN * * * 

* TOLERANCE * +** 



*************** 

* * 
♦COMPUTE FIRST* 

* TERM OF * 

* EO.UI * 

* * 
*************** 



* COMPUTE NEXT* 

* TERM OF * 

* EQ.(2» * 

* * 
*************** 



*************** 



*************** 



************** « 



* *YES 

IS DIFFERENCE * 
BETWEEN TFRNS * 
LESS THAN * 

* TOLERANCE * 



•YES **• 

* * * 
K * 30 •) CI* 



********** )* go* 



*******:***)* CO* 
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SUBROUTINE BESK 



*************** 
ENTER 

*************** 



*************** 



*HO***>* COMPUTE KO 
* * * KO+TERM 



*************** 



*************** 



*CO***)* COMPUTE KQ. * 

* * *K1 FU.(2),<3)* 

** * * 

*************** 



♦ NO 
* 
THAN 



********* 



♦YES ** 

* * * 

TERMS *********#ci* > 



********* 



10 v 
** ************* 

* s;.t fsuok * 

* Ci..n(= v:i i * 
^ * 
*************** 



*************** 

* * *** 

* COMPUTE NEXT* * * 

* TERM OF ****) 80* 

* EQ.(7> * * * 

* * *** 
*************** 



*************** 

* * 

* SET ANSitfEK TO* 

* BO OK ,11 FOR* 

* N = CJK 1 * 

* * 
*************** 



**************** ***. 



*************** 

******** v****** 



*************** 

* SUBTRACT * 

* COMPUTED * 

* CONSTANT TU * 

* COMPLETE * 

* EQ.m ♦ 
*************** 



*************** 

rxiT 

*************** 



****#************** 



***************** 



********* 



*************** 

* * 

♦COMPUTE FIRST* 

* TERM OF * 

* E0.C9). SET * 

* Kl = * 
*************** 



*************** 



# K0| K R = M 

* 
************** 



* * ** *********** 



*************** 



*************# 



S ********* 3> 



**rt>i***>* COMPUTE Kl «■* 
** * * Kl+TERM * 



*************** 



*C5***)+ CUMPUTt" TEnM* 

* * * L OF t-i.i II * 

** * * 

*************** 



********** 
FXIT 

****>**** ** 



♦ YES 

ft TERMS 



********* 



********* 



******************* 



*************** 

■■? # 

* SFT purou * 

* C'.iOE TO O (NO* 

* FftiUJR ) w 

* * 
******** ****** v 



*************** 

* * *** 

* COMPUTE NfXT* * * 

* TEftM OF #***) d5* 

* FQ.I9I * * * 

* * *** 
*************** 



*************** 

* * 

* * 

* STEP L 



*************** 



* *** 

* * * 
****) Ci>* 

* * * 
*** 



******************* 



******************* 



******** 



*************** 

* * *** 

* ADD CQMPUTEO* * * 

* CONSTANT TO ****> CI* 

* CQMPLETF * * * 

* EQ.19) * *** 
*************** 



*************** 

* * 
♦SET ANSWER TO* 

* LAST VALUE * 

* FUfc E(J.( 1) * 

* * 
*************** 



***************** 



"16 V 

************ **# 

* * *** 
♦COMPUTE FIRST* * * 

* TERM OF ****) BO* 
+ F0.( 7>. SFT * * * 

* KO = + *** 
*************** 



25 V 

*************** 

* * *** 

* COMPUTE GO * * * 

* ANO 61 FOR X****l CO* 

* Ea.14).<5> * * * 

* * *** 
*************** 

105 



*************** 

EXIT 

*************** 



SUBROUTINE CILI1 



*************** 

ENTER 
*************** 



*************** 



* 15 -CK LESS * *RESULT EQUALS* 
THAN OR EOUAL ***************)* INFINITY. * 

* TO 1.0 * * IER-1 * 



*************** 



* * *************** 

* INITIALIZE * * 

* ARITHMETIC * * EXIT 
•AND GEOMETRIC* * 

♦MEANS ARI.GEO* *************** 

*************** 



* * * SET AARI^ARI* 
«.A4***I* COMPUTE * 

* * * ARI*ARI+GEO * 
** * * 

*************** 



* * *************** 

*YES * • 

* * * 

IS ARI-GEO ***************)*RESULT EQUALS* 

NEGLIGIBLE * * PI/ARI * 

* * * 

* * *************** 



* COMPUTE NEW * 

* GEOMETRIC * * EXIT 

* MEAN GEO * * 

* * *************** 
*************** 



V 30 
*************** 

* * *** 

* * * * 

* REPLACE ARI ****] (U* 

* 8Y ARI/2 * * * 

* * *** 
*************** 



SUBROUTINE CELI2 



*************** 

********** ***** ** * * ** * * 

• * * * SET AA(U«ARI* * * * * 

FNTFR * *(vo***)* COMPUTE * *CO*****RESULT EQUALS************ 

^' fcK * * * * ARI-ARI+GEO * * * * -INFINITY * * 

*************** ** * * ** t^^^^^^^^^l I 

* *************** *************** * 



*************** 

* ** * * 

* * * * * 
I IPR=0 * * COMPUTE * *C1*** »*RESULT EQUALS* 

* * * AN=W/ARI+AN * * * * ft * 

* * * * ** + * 
*************** *************** *************** 



*( ***************** 



* *YES * *YES 
i <; rtc t f<is * * * * 

THAN OR EQUAL ********* * IS ARI-GEO ********* * EXIT 

* ™ LO * * * # NEGLIGIBLE^ * * 



*************** * 

* * * * ¥ *»<* 

* * * * * * * 

* IFR = l # * * COMPUTE NEW ****) A9* 

* * . * * GEO * * * 

* * * * * *** 
*************** * *************** 



*************** 



IS B LESS 
THAN 



******************* 



*************** * * 

* * * *YES *** 



IS fl EQUAL TO 




*************** 



* * 10 


30 V 


* * 


*************** 


♦YES *** 


* * 


* * * 


* * 


[S GEO EQUAL *» B4* 


♦RESULT EQUALS* 


TO 0.0 * * * 


* INFINITY * 


* *** 


* * 


* * 


******* ******** 



* INITIALIZE * * _„ 

* ARITHMETIC * * EXIT 

* MEAN* ARI * * 

* n "™» AKl * *************** 

*************** 



******************* 



*************** 

* * * » 

* * * * 

* INITIALIZE * *RESULT EQUALS* 
*AA, AnJ AND W* * PIM*AN/ARI * 
*************** *************** 



60 V 

*************** 
** * * *** 

* * * COMPUTE * * * * 
j,At)***|* W= 2(H+AA*GE0>+**+) 80* * EXIT 

* * * SET AA^AN * * * * 

** * * *** *************** 

*************** 



SUBROUTINE EX PI 



*************** ** * *YES 

* « * * J 

ENTER * *BO*l IS X EQUAL TO ********* 

* * * * o * * 
*************** ** * * * 



I ... „ * * EVALUATE * 

* [ER»0 * * POLYNOMIAL ♦ 

* * 'APPROXIMATION* 
*************** ******.****;»*»J 



.♦,.,„*. : RE ^tyNoSYa s : 

* * VALUE PLUS * 

* *-ALOG(ABS(X))* 

* *************** 



♦NO *** * *************** 

* * * * * 

IS X LESS *) BO* * * EXIT 

THAN -«, * * * * * 

* *** • *************** 



******************* 



* IFR-l * * * RESULT IS * 

J J INFINITY * 

*************** « **********.***J 



EXIT * * * EXIT 

* * * 
*************** * ********»****«» 



******************* 



* SET ARGUMENT* 
♦EQUAL TO 4./X* 

* * 
*************** 



*************** 

* * 

* EVALUATE * 

* POLYNOMIAL * 
♦APPROXIMATION* 

* * 
*************** 



*************** 

* * 
♦RESULT EQUALS* 

* POLYNOMIAL * 

* VALUE TIHFS * 

* EXPI-XI * 
************ *** 



EXIT 

*************** 



SUBROUTINE SIC1 



*************** 



*************** 



*************** 



* SET Z EQUAL * 

* TO ABS(K) * 



*************** 



1 GREATER 
THAN 4 



*************** 



*************)* REPLACE Z BY* 
* 2*2 * 



*************** 



******** ******* 
w * 

* SET SI EQUAL* 
*TO SINIZ) AND* 

* Y EQUAL TO * 

* COS til * 
*************** 



*************** 

* . * 

* EVALUATE * 

* POLYNOMIAL * 

* EXPRESSION * 

* FOR SI * 
*************** 



*************** 



* REPLACE 1 BY* 

* 4./Z * 



*************** 



IS ARGUMENT 
EQUAL ra 



***************}* ci EQUALS 
* INFINITY 



****** ********* 



*************** 

* EVALUATE * 

* POLYNOMIAL * 

* EXPRESSIONS * 

* FOR U AND V * 
*************** 



* EVALUATE * 
♦APPROXIMATION* 

* FOR CI * 

* * 
*************** 



*************** 



*************** 



*************** 

* * 

* CALCULATE * 

* CI=Z*<Sl*V * 

* +Y*U) * 

* * 
*************** 



*************** 



*************** 



*************** 

* * 

* CALCULATE * 

* SI— Z*ISI*U * 

* +Y*V) * 

* * 
*************** 



*************** 

* * 
+AOJUST RESULT* 

* VALUES FOR * 

* NEGATIVE * 

* ARGUMENT X * 
*************** 



*************** 



*************** 



SUBROUTINE CS. 



*************** 

ENTER 
*******#*#*++** 



*************** 



* SET Z EQUAL * 

* TO ABSIX) * 

* * 

*************** 



* *************** 
*YES * » 

* * SET D EQUAL * 

Z GREATER ****#**********)*to COSIZI AND* 

THAN * • * S EQUAL TO * 

• * SIN(Z) * 

* *************** 



*************** 



EVALUATE 
POLYNOM" " 



* CALCULATE * * PSLYNOHIAL » 

* C*SORTIZ> * * EXPRESSIONS * 

* * * FOR A ANO B * 
*************** ************«*, 



* SET S EQUAL * * * 

* TO Z*C * * REPLACE I BY* 

* REPLACE Z BY* * SORTIZI * 

* Z*Z * * * 
*************** *************** 



*************** 

* * * * 
*. EVALUATE * * CALCULATE * 
♦APPROXIMATION* * C-.5* » 
•FOR C ANO FOR* * Z*( D*A+S*B> • 

* S * * * 
*************** *************** 



«„,, * * CALCULATE * 

EXIT * * S-.5+ * 

* * Z*(S*A-0*BI • 

*************** * * 

*************** 



EXIT 

*************** 
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SUBROUTINE RTMIT 



**#*#********** ** *************** 



*************** ** *************** 



IS ABS(FX-X) * 
LESS THAN ********* 



* NO * * 

**********)* 36* 
* * 
*+* 
11 



*************** * * 

** * *YES 



l™?9.... * *B1***)* ITERATES * *Cl*> IS I»ITEND ********* 

m&ve^ 



* A=FUNCTIY) * * * * VALUE ANO 

* * ** * HOVE IT IN A* 
*************** *************** 



* *NO * *YES * * *** ** 
DOES Y SOLVE * * * * ***** 
THE EUUATION ********* * IS X NEARLY ********* * ADD 1 TO I ****) A6» ** 

GIVEN PY * * * EQUAL TO * * * ***** 

* FUNCT * * * * * * * *** ** 
* + * * * * *************** ** 



******************* 



* * * * *YES *** * * * 

* * * * IS RELATIVE * * * * * * 

1* X=A * * * ERROR LESS *) C4* * * IER = l * 

* * * * THAN EPS * * * * * * 
+ * * * * *** * + * 
*************** * * * 20 * *************** 

* * * * *** * * 

* * * NO * * * * 

* * **********)* ex* * * 



*( ***************** 



V * *************** 

*************** * * * *1 

* * * * * * * 

EXIT * * * COMPUTE * *C4*)* EXIT 

* * * FO-FUNCTIOI * * * * 
*************** * * * ** *************** 

* *************** 



******************* 



* * * DOES X=0 * 

* SET 1=1 FOR * * SOLVE THE ********* 

* ITERATION * * EQUATION * * 

* * * * * 
*************** * * * 



* COMPUTE * * * * * 
>* ITERATED * *B6***|* X-0 * 

* VALUE A * * * * * 

* * *# * * 
*************** *************** 



*N0 *** 

IS NEH A * * * * 

DIFFERENT *» A3* * EXIT 
FROM OLD A * * * * 

* *** *************** 



********** **** ***** 



* DOES *ND *** 
ITERATION * * * 

MOVE IN AREA *» Bl* 
DERIVATIVE * * * 

* NEARLY 1 * *** * * 
* * *************** 



* x=A IFR=2 ** 

* * 

* * 
*************** 





♦ YES 


*** 


IS 


ABSIFXI * 


* * 


GREATER THAN 


*l CI* 




EPS * 


* * 




*■■ 


*** 


* 


* 




: 


* * 


*** 




* NO 


* * 




**********}* CO* 






* * 






*** 






14 



SUBROUTINE HTHIT 



*************** 

ENTER 
*************** 



*************4i* 
** * * 

*. * * PERFORM * 
*BD***|*ITERATION FOR* 
* * * H»ltITEND*l * 

** * * 

*************<H* 



*VES ' 

* * * * 

K* I TEND ****+*****DO*** 1 * 



****** ********* 



*************** 



*************** 



IER*0 * 

FL=FUNCT<XLI* 



*************** 



*************** 

* * 

* PERFORM * 

* BISECTION FOP:* 

* K-l.lTEND * 

* * 
*************** 



*************** 

* * *** 
+ * * * 

* ADD 1 TO K ****) 62* 

* * * * 

* * *** 
*************** 



******************* 



* *N0 
DOES XL SOLVE * 

THE EQUATION 
RIVEN BY * 

* FUNCT * 



*************** 

** * * 

*********I B2 t**|t B ^S^^ON E NiJJ* 
** * * x-VALUE AND * 

* ** * F-FUNCTI-XI * 

* *************** 



*************** 

* * 

* * 

* SET IER-2 * 

* * 

* * 
*************** 



* *VES 

IS F*FL * 
GREATER THAN 



********* 



*************** 



* X = XL 

* * 

* + 

*************** 



*** 

* * * 
**•*) D6* 

* * * 

*** 



DOES X SOLVE 
THE EQUATION 



♦YES *** 

* * * 

*i 06 : 



********* 



************** * 



*************** 



*** 

* * 

D5* 



******************* 



******************* 



*************** 

* * 

* FR»FUNCT(XR»* 

* * 

* * 
*************** 



*************** 



IS F*FR LESS 
THAN 



********* 



* SET X=XL 

* * 

* * 
*************** 



* *** 

* * * 
***♦) D6* 

* * * 
*♦* 



*************** 

* * 

* SET XL*X , * 
*FL=F t XK=XM»* 

* FK=FM * 

* * 
*************** 



************** *#*** 



DOES X« SOLVE 
THE EQUATION 



********* 



*************** 

* * 

* INTERCHANGE * 

* XL WITH XR * 

* AND FL WITH i> 

* FR * 
*************** 



*D5***)* ADO 



*************** 



1 TO N #***) H2* 

* * # 

+ ** 



*************** 



*( ***************** 



******************* 



*************** 



*************** 



91 

#** 

* * * 
****) 06* 

* * * 
*** 



IS ABS(f) 

LESS THAN 
AbS(FR) 



* * *CDMPUTE OX BY* 

*******+** C6 *+*]* INVERSE * 
QUADRATIC 



91 * 

V 

*************** 

EXIT 

*************** 



******************* 



*VFS *** 
* * * 
*) BO* 

* * 
*** 



* * 3i 


* 


V 


* * 


* 


*************** 


1 *YES *** 


* 


* + 


* * * 


* 


* * 


IS N LESS ♦) C6* 


* 


* SET XM=X AND* 


THAN ITENO+1 * * * 


* 


* FM=F * 


* *** 


* 


* * 


* * 


* 


*************** 



*( ***************** 



*************** 



* CALCULATE * 

* X-{XL+XRI/2 * 

* * 
*************** 



SET XR-X AND* 
FR*F * 



*************** 



*************** 

* * 

* * 
♦COMPUTE NEW X* 

* * 

* * 
*************** 



*************** 



*** 

* * * 
IER»1 ****> 06* 

* * * 
*** 



*************** 



21 

S *** 

* * 
*> C3* 

* * 
*** 

* 
* *** 
* NO * * 
**********)* CO* 



BISECTION 
STOP 



112 



91 

♦YES *** 

* * * 
*> 06* 

* * * 

* *** 

*** 

* * 
********** }* oo* 

* * 
*** 

34 



it£rXt?on 

STOP 



NO 



SUBROUTINE RTNIT 



*************** ** * *YES 

» * * * COULD * 

FNTFR * «BO») ITERATION ********* 

* * * * STOP * * 

*************** ** * * * 



* N=Ot IFR-O, * * * * * * 

* PP.EP4RE * *B1*I IS N tESS^ »***»*♦*» 

* ITFRATION * * * ' ~ " - 

* LOOP * ** 
******** ******* 



V V 

*************** *************** 

** * * * * *** 

* * * » PREPARE FOR * * * 

«A>***|* »DO 1 TO N * * NEXT ****) A2* 

« * * * * ITERATION * * * 

** * * * * *** 

*************** *************** 



******************* 



«. *NO * * 

DOES Y SOLVE * * NO * 

THE FQU4TION ********* * CONVERGENCE.* 

GIVEN BY * * * IER=1 * 



FUNCT 



*************** 



*( ***************** 



1 io V * 150 * 

*************** * V 

* * * ** *************** 



* * * ** *************** 

*************** * 



*************** * ** * *N0 

* ** * * IS X=0 ROOT * 

FXIT * **B5*I OF THE ********* 

* «* * * EOUATION * * 
*************** * ** * * * 



******************* 



IS DERIVATIVE ********* 
EOUAL * * 



*************** 



120 V 110 * * 

*************** * V 

* * *** * *************** 

* ITERATION IS* * * * * , „ 
*NOT POSSIBLE.****) A«* * * EXIT 

* SET IEK=? * + + * * 

* * *** * *************** 
*************** * 



******************* ******************* 

* * 

* V 
70 * + 150 

* .* 

* * * *¥ES *** 
♦COMPUTE NEW X* * COULD * * * 

* BY NEWTONS * * ITERATION *l B*J 

* HETHnD * * STOP * * * 

* * * * *** 
*************** * * 

* * * *** 

* * NO * * 

* **********)• Bl* 







* * 




50 




# 


* 






* 






*YES 

* 


*** 
* * 




rs 


ABS(X) 




*) 85* 


LESS 


THAN EPS * 


* * 


* 


* 


* 
* NO 


* 


*** 

*** 
* * 






**********)* BO* 










* * 










*** 










80 



SUBROUTINE PCILRT 



****«<********** 
ENTER 

*************** 



*•*•«*»•******* 
** • * 

* * * SOLVE 111), • 
•BO«**l* 1121 FOR * 

* * * CORRECTION * 
** * TERMS * 

*************** 



THE EFFECT OF 

THE LOGIC IS 

TO REFINE 

ROUT BY 



SET ERROR #****♦****** 
CODE TO 4 » * 



*************** 



*************** 



* SET ERROR * 

* CODE TO * 



*************** 



******* ****** ** 
♦FORM MODIFIED* 

* ROOT BY * 

* APPLYING * 

* CORRECTION * 

* TERMS * 
*************** 



DOING FINAL 

ITERATIONS ON 

ORIGINAL 

COEFFICIENTS 



*************** 



*D1***>* SET ERROR ********* 
* * * CODE TO 1 * 



*************** 



* *YES *** 

* IS HIGH ORDER * * * 

COEFFICIENT *l DO* 

* ZERO * * * 

* * *+* 



ARE * * * 

CORRECTION **********C2*| 
TERMS LESS * ** « 
♦THAN .OOOOL* * ** 



120 * * 

*IS RATIO OF*YES 

* COMPLEX TO 
REAL ROOT 

* LESS THAN 
• .0001 



**********02**« ) * 



SET ERROR 
CODE TO 2 



*************** 



10 * * 15 

* * 

* *YES *** 

IS OROER LESS *1 Dl* 
THAN 1 * * * 

* * *** 



*************** 

* * 

* STEP * 

* ITERATION * 

* COUNTER * 

* * 
************* ** 



*(*••** 4* ********** 



*************** * 

* * * : 

* PREPARE FOR * ** 

* SYNTHETIC *******#*****03*)« 



*************** 



*************** 



*************** 



******************* 



32 * 



30 



* *YFS *** 
IS ORDER * * * 

GREATER THAN *I 02* 
36 * * * 

* * *** 



* * 60 
* * 

*N0 *** 

HAVE 500 * * * 

ITERATIONS *l A9* 

BEEN DONE * * * 

w *** 



♦FORCE ROOT TO* 
* BE REAL * 



*************** 



*( **************** 



40 V 
*************** 
♦MOVE ORIGINAL* 

* COEFF jn * 

* REDUCED * 

* COEFFICIENTS* 

* ARRAY * 
*************** 



90 



50 



* HAVE FIVE <-N0 *** 

* PAIRS OF * * * 
INITIAL •) A7* 

* VALUES BEEN * * * 

* TRIED * *** 



.,». 



* PERFORM * 

* SYNTHETIC * 

* DIVISION * 

* * 
*************** 



*************** 

* * 

* ASSUME * 
I* INITIAL * 

* VALUES FOR * 

* ROOTS * 
*************** 



* SET ERROR **•*) 03* 

* CODE TO 3 * * * 

* * *** 
*************** 



*************** 

* FORM NEW * 

* REDUCED * 

* COEFFICIENT * 



******************* 



** 



*************** 



* * * INCREMENT * 
*A7***I* INITIAL * 

* * * VALUES AND * 
** * COUNTER * 

*************** 



* SWITCH TERMS* 

* IN REDUCED * 

* COEFF ARRAY * 
♦WITH ORIGINAL* 
**•♦♦♦********* 



HAS ROOT 
COMPLEX 



********* 



*************** 
** * * 

♦ * * INITIALIZE * 
*A8***I* TEMPORARY * 

* * * VALUES * 
** * * 

*************** 



* IS THIS * * • 

* ♦ coePPYcient * ** A8 J 

* ARRAY * *** 

* * 

* * *** 

* YES * * 

*******«**)* C 2* 



* FORM COMPLEX* 
♦CONJUGATE AND* 
♦PlACFlN ROOT* 

* TABLE • 
****** ********* 



* I ***************** 



60 V 

*************** 
** * SOLVE * *** 

* * * (31,1*1, * * * 
*A9***)* (7)-(l0) FOR****» BO* 

• * * NEWTON * * * 
** *RAPHSON TERMS* *** 

*************** 



165 * * 45 

• * 

* *YES *** 

* * * * 
♦ ARE ANY ROOTS *) A6* 

* REMAINING * * * 

* * *** 

* • 

* * *•* 

* NO * * 

***»******)« 03* 



SUBROUTINE PAOD 



*************** 

FNTFR 
*************** 



*************** 

* * 

* DEFINF * 

* DIMENSION OF* 

* RESULTING * 

* POLYNOMIAL * 
*************** 



*************** 

* * 

* -AOD * 
♦CORRESPONDING* 

* COEFFICIENTS* 

* * 
*************** 



*************** 



*************** 
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SUBROUTINE PACDM 



*************** 

ENTER 
*************** 



* DEFINE * 

* DIMENSION DF* 

* RESULTING * 

* POLYNOMIAL * 

*************** 



* COMBINE * 
•CORRESPONDING* 

* COEFFICIENTS* 

* * 

*************** 



*************** 

EXIT 
*************** 



SUBROUTINE PCLA 



*************** 
ENTFR 

** ************* 



*************** 



+ SET * 

* IDIMY-IOIMX * 



*************** 



*************** 

* * 

* Move * 

* POLYNOMIAL X* 

* INTO * 

* POLYNOMIAL Y* 

*************** 



*************** 



***********:**** 



11 V 



SUBROUTINE PSUi 



*************** 



*************** 



* DEFINE * 

* DIMENSION OF* 

* RESULTING * 

* POLYNOMIAL * 

*************** 



* SUBTRACT * 

•CURRESPONOING* 

* COEFFICIENTS* 

* * 
*************** 



*************** 



*************** 
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SUBROUTINE PHPY 



*************** 

ENTFR 
*************** 



*************** 

* * 

* DEFINE IOIMZ* 

* AS DIMENSION* 

* OF RESULTING* 

* POLYNOMIAL I* 
*************** 



* 400 PRODUCTS* 

* xi Li*rui to* 

* FORM ZIK) * 

* * 
*************** 



*************** 



SUBRQUTIHE PC'IV 



*************** * 

* * * * 

ENTER * »bo**»)* SET ERROR 

* * * * CDOE TO 

*************** ** * * 

* *************** 



V 40 * 

*************** V 

* PNtlRM * ** *************** 
******* ******* * * * 

* NORMALIZE * *R1*1* EXIT 

* POLYNOMIAL * * * * 

* * ** *************** 
*************** 



* *NQ 

DIMENSION OF * 
V LARGER THAN ********* 

ZF.'<0 * * 



* *YES * 
I3IM6NSION OF * * 

V LARGER THAN ********* 

DIMENSION OF * ** 

* X * *« 
* * ** 

* * *+ 



*************** 



* xt ii/yi n * 



*************** 



*************** 

* * 

* SUBTRACT * 

* MULTIPLE OF * 

* DIVISOR * 

* * 
*************** 



*************** 



*************** 



+ ** 

* ** 

90 V 30** 

*************** ** 

♦PNORM * *** ** 

******* ******* * * ** 

* NORMALIZE ****! BO* ** 

* REMAINDER * * * ** 

* * *** ** 
*************** ** 

** 
** 



******************* 

* 

* * 

50 V 40* 

*************** * 

* * *** « 

* * * * * 

* SET ERROR *) &1* * 

* cooe TO I * * * * 

* * *** * 
*************** * 



****************** 



*************** 



*************** 



SUBROUTINE PASO 



*************** 

ENTER 

*************** 



*************** 

* * 

* INITIALIZE * 

* A-B«Oi * 

* JalDIMX * 

* * 
*************** 



* * 3 

* 

*NO *************** 

* * 

J GREATER *************)* EXIT 
THAN * * 

* *************** 



* SET Z=P*A+B * 

* AND * 

* B=Q*A+X(JI * 

* * 
*************** 



* REPLACE 

* Z 



*************** 



* SURTRACT 1 ****! A2» 

* FRUM J * * * 

* * * ** 
*************** 



SUBROUTINE PVM. 



*************** 



*************** 



*************** 



********* 



* PERFORM * 

* NESTED * 

* MULTIPLI- * 

* CATION * 
*************** 



*(***************** 



1 * 
V 
*************** 

EXIT 

*************** 



SUBROUTINE PVSUB 



Or************** ** * * 

* * * * 

ENTER * *BO***>*SET 0IHENSION 

* * * * OF I TO ZERO 

*************** ** * * 

* *************** 



* * I 2 * 

* * y 

* *Y£S *** ** *************** 

* * * * * * 

IS X A ZERO *) BO* *Bl*)+ EXIT 

POLYNOMIAL ****** 

* * *** ** *************** 



* * 2 

* * 

* *YES *** 

* * * 
IS Y A ZERO *) Bl* 

POLYNOMIAL * * * 

{• * *** 



*************** 



V- 
*************** 

** *PHPY * 

* * ******* ******* 

*A4***)*HULTIPLY Y BY* 

* + * 1 * 
** # * 

*************** 



V 
*************** 
*PCLA * 

******* ******* 

* MOVE WORK * 

* AREAS * 

* * 
*************** 



* SET FACTOR = 

* XUDIMX-I) 



*************** 



*************** 

♦PADOM * 

******* ******* 

* ADO I TO * 

* FACTOR*WORK * 

* * 

*************** 



♦YES *** 
* * * 
*> Bl* 



* ADD 1 TO I ****) A4* 

* * * * 123 

* * *** 

*************** 



SUBROUTINE PCi.0 



*************** 

ENTFR 
***** ********** 



* SET K-l * 

* * 

* * 
*************** 



*A2***>* SET J-IOIMX * 

* * * * 

** * * 

*************** 



***************)* A [)0 1 TO K 

* 
************** 



* * * *YES *** 

* REPLACE * * IS IDtMX * * * 



** * ************ 



V 2 
*************** 




5 * 
V 


* * *** 




*************** 


* * * * 


# 




* SUBTRACT 1 ****) A3* 


* 


EXIT 


* FROM J * * * 


* 




* * *** 




************** * 


*************** 







SUBROUTINE PILO 



*************** 



*************** 



* INITIALIZE * 

* OUADRATIC * 

* SYNTHETIC * 

* DIVISION * 

*************** 



COMPUTE VALUE 

OF DERIVATIVE 

BY SUBROUTINE 

t POSO * 

*************** 



*************** 



*************** 



EXIT 

*************** 



SUBROUTINE POEft 



*************** 

ENTER 

*************** 



* *NO * * 
IS IOIMX * * RESULT IS * 

GREATER THAN ***************J* ZERO * 

1 * * POLYNONIAL * 

* * * * 
* * *************** 



*************** 



*************** 

* * 
+ * 

* SFT I=*D * 

* * 

* * 
*************** 



EXIT 
*************** 



* ft «V***)* ADO I TO I * 

* * * * 

** * * 

******** ******* 



*************** 



*NO *** 

* * * 

I EQUAL TO *l A4+ 

IDIMY * * * 

* *** 



*************** 



SUBROUTINE PINT 



*************** 



*************** 



*************** 

* * 

* SET * 
*IOlMY»If>IMX+i* 

* AND YUI-0 * 

* * 
*************** 



GREATER THAN ********* 

* * 



*************** 



«M***>* ADO I TO I * 

-. * * * 

** * * 

*************** 



* SET * 

♦YUtU-XIII/t* 



*************** 



is i equal to 
ioihx 



* t ***************** 

* 
* 
1 



V 

*************** 



*************** 
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SUBROUTINE PfiEO 



*************** 
ENTER 

*************** 



I V 

*************** 

*# *PDIV * 

* ******* ******* 

» ai***)*d[vioe x ev y* 

* * * * 

** * * 

** ****** **** *** 



********* 



* *NQ * 

DIMENSION OF * * 

X LARGER THAN ********) 

ZERO * * 



*************** 



V 1 

*************** 

* + *** 

* INTERCHANGE * * * 

* DIMENSIONS X****» At* 

* AND Y * * * 

* * *** 
*************** 



******************* 



*************** 
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SUBROUTINE mOMt 



ENTER 



I * * 

* * 

• NO 



gre»ter'th*n " *«*«**»***»*»*********** 
o * I 



« * v 

, » N Q *************** 

IS XCIDIHX) *»*****«»****»* EXIT 
■ NEGLIStOLE * * ,„„„„.»„„ 



1 

*•* 

• * * * 
♦REPLACE [OIHX«***l »1« 
« BY IDIMX-l * * * 

* * **♦ 
*************** 



